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Abstract 

We consider a hybrid diffusion process that is a combination of two Ornstcin-Uhlcnbcck 
processes with different restraining forces. This process serves as the heavy-traffic approximation 
to the Markovian many-server queue with abandonments in the critical Halfin-Whitt regime. We 
obtain an expression for the Laplace transform of the time-dependent probability distribution, 
from which the spectral gap is explicitly characterized. The spectral gap gives the exponential 
rate of convergence to equilibrium. We further give various asymptotic results for the spectral 
gap, in the limits of small and large abandonment effects. It turns out that convergence to 
equilibrium becomes extremely slow for overloaded systems with small abandonment effects. 
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1 Introduction 

Within the fields of stochastic processes and queueing theory, the Halfin-Whitt regime refers to 
a mathematical way of establishing economies-of-scale in many-server queueing systems like call 
centers (see [13]). The Halfin-Whitt regime in fact prescribes a scaling under which the many- 
server systems converge to limiting processes, which are for most systems diffusion processes. This 
paper deals with many-server systems in the Halfin-Whitt regime with the additional feature that 
customers are impatient, and may abandon the system without being served. For such systems with 
abandonments, we are interested in the spectral gap, which is inversely related to the relaxation 
time or the speed at which a system reaches stationarity. A large relaxation time in general indicates 
that replacing time-dependent characteristics by their stationary counterparts might lead to poor 
approximations. As it turns out, the rate at which customers renege (abandon the system) greatly 
influences the relaxation time. 

In recent years, a large number of papers have dealt with the influence of reneging or abandon- 
ments on the system behavior (see, e.g., [8, 15, 35, 36, 37, 41, 42]), and it is widely accepted that 
reneging is indeed one of the main factors driving the system performance. One of the key insights 
is that the system behavior strongly depends on whether it is stable or overloaded. By stable we 
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mean that it can serve all customers, even if none of the customers would abandon the system. The 
spectral gap (or relaxation time) is also very different for stable or overloaded systems. In fact, 
we find that stable systems have a relatively short relaxation time, whereas the relaxation time of 
overloaded systems can become extremely large, particularly when the reneging rate is small. 

The model we shall consider is the M/M/s + M system, better known as the Erlang A model. 
This model is a standard Markovian many-server queueing system with Poisson arrivals, exponential 
service times, s servers, and with the additional feature that customers that are waiting in the queue 
abandon the system after exponentially distributed reneging times. The queue length process in 
the Erlang A model, denoted by (Q(t))t>o, is a birth-death process. Whitt [36] (see also [38] and 
[19]) derived a fluid approximation for the the steady-state behavior of the overloaded Erlang A 
model, and he further showed that a diffusion limit might provide refined approximations. Garnett, 
Mandelbaum and Reiman [15] proved a diffusion limit for the Erlang A model in the critical regime. 
In particular, they showed that under certain conditions a sequence of normalized queue length 
processes converges to a certain diffusion process (X(t))t>o- These conditions are in fact the ones 
that correspond to the Halfin- Whitt regime, in which the arrival rate A and the numbers of servers 
s are scaled such that, while both A and s increase toward infinity, the traffic intensity po = X/s 
approaches one, with 



The diffusion process (X(t))t>o is a combination of two Ornstein-Uhlenbeck (OU) processes with 
different restraining forces, depending on whether the process is below or above zero. The number 
of customers in the Erlang A model can be roughly expressed as s + yfsX(t) for s sufficiently large. 
The diffusion process is generally easier to study than the birth-death process, and can thus be 
employed to obtain simple approximations for the system behavior. The steady-state distribution 
of the diffusion can be easily obtained (see (6) below), but less is known about the time-dependent 
behavior. In this paper we shall present an explicit and asymptotic characterization of the spectral 
gap of (X(t))t>o- The spectral gap of the diffusion process provides an understanding of the 
relaxation times for the Erlang A model in the Halfin- Whitt regime. Most importantly, we shall 
study in detail the impact on the spectral gap of the capacity parameter f3 and the reneging rate rj, 
which shall enhance our understanding of how the Erlang A model behaves for positive/negative /3 
and small/large values of rj. 

The diffusion process (X(t))t>o also applies to the G/M/s + M system, in the same asymptotic 
limit, which was proven by Whitt [37]. Stochastic processes for more general systems with aban- 
donments were obtained recently by Dai, He and Tezcan [8] for the G/Ph/s + M system. In this 
case, the limiting process is still a diffusion process, but it becomes multi-dimensional. Zeltyn and 
Mandelbaum [42] derived approximations for the M/M/n + G in the Halfin- Whitt regime. In case 
of general service times, the limiting process is not even a diffusion process (see e.g. [27, 37] for cases 
without reneging). Therefore, the one-dimensional diffusion process {X(t))t>o strikes the proper 
balance between simplicity and tractability, while retaining the essential features of abandoning 
customers in many-server systems. 

The Erlang A model is particularly interesting, as it incorporates three classical queueing sys- 
tems as special cases. In the case of no reneging (with rj = 0) the Erlang A model reduces to the 
Erlang C model, or M/M/s system. Halfin and Whitt [17] established that the limiting process 
behaves as a Brownian motion above zero and an OU process below zero. In [23] we have referred 
to this process as the Halfin- Whitt diffusion. For rj = 1 the Erlang A model becomes an infinite 
server queue or M/M/oo system, for which the stochastic-process limit is known to be an OU 
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process [18]. For rj — > oo the Erlang A model becomes the Erlang B model or M/M/s/s system, in 
which case the stochastic-process limit is a reflected OU process (see [12, 24, 35]). For the diffusion 
approximations, we show the rj — > oo reduction in Section 3.3. Our analysis of the spectral gap of 
(X(t))t>o provides results for each of these three cases. 

Mathematically, determining the transient distribution for the present diffusion process involves 
analyzing a Schrodinger type equation (see Section 4) with a piecewise parabolic potential function, 
or, equivalently, a Fokker-Planck equation with a piecewise linear drift (see (2) and (3) below). Such 
problems arise in a variety of other applications, such as linear systems driven by white noise [3, 2], 
the Kramers' problem [26] and escape over potential barriers [22]. Invariably, the solution involves 
the parabolic cylinder functions, and these we discuss in detail in Section 5. 

The key in determining the spectral gap is in fact determining the Laplace transform of the 
transient probability distribution over time. The spectral gap then follows from the dominant 
singularity of the Laplace transform. The main results are presented in Section 2, the three special 
cases (rj = 0,1, oo) are discussed in Section 3 and the proofs are given in Section 6. Before the 
proofs we give some basic background on Schrodinger equations (Section 4) and parabolic cylinder 
functions (Section 5), whose properties are heavily used later. In Section 7 we establish monotonicity 
properties of the spectral gap. 



2 Main results 

The diffusion process (X(t))t>o is a Markov process on the real line with continuous paths and 
density p = p(x, t) = p(x, t; xq] P, rf) that satisfies the forward Kolmogorov equation 



where 



:?:r III (3) 



and (with 5(-) the Dirac function and p x = dp/dx) 

p(x,0) = 5{x-x ), (4) 
p(0+,t)=p(0-,t), p x (0+,t)=p x (0-,t), (5) 

and p(x,t) must decay as x — > ±oo. The limiting distribution of the diffusion process is (see [15]) 

6 -ix 2 6 Br ' X n ( 6 ) 

e a x e px , x < 0, 

where C~ l = J °° e'^dx + e^V^dx. 

As shall be discussed in Section 4, this problem has a purely discrete spectrum for all r\ > 0, 
and it is confined to the real axis. The spectral gap can thus be defined as the absolute value of the 
least negative eigenvalue of the operator in the right-hand side of (2). It governs the asymptotic 
rate of convergence to the stationary distribution. An alternative description of the spectral gap 
is the absolute value of the singularity closest to the imaginary axis in the range Re(0) < of 
the Laplace transform p. Denote this dominant singularity by and the spectral gap by r(/3,r]). 
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The relaxation time, which measures the time it takes for the system to approach its steady-state 
behavior, is defined as (see [4, 7]) 

r = inf{T : p(x, t; x ; /5, rj) - p(x, oo; x ; (3, r)) = 0(e~ t/T )}, (7) 

and hence r _1 = — Re(#) = r(fi,rj). For this problem 9 is real, so that —9 = r(/3, 77). Our definition 
of the relaxation time in (7) assumes the initial condition p(x,0) = 5(x — xq) in (4), and then the 
approach to equilibrium is governed by Ai = r. But we could certainly have initial conditions that 
would lead to a faster approach. For example, if p(x, 0) = p(x, 00) then p(x, t) = p(x, 00) for all 
t and equilibrium is attained instantaneously. We could also have initial distributions p(x, 0) that 
have zero projections on, say, the first L eigenfunctions, and then the sums in (31) and (32) below 
would be replaced by Y^=L+i e ~ Xntc n^n{x), where the c n may be computed in terms of p(x,0). 
Then the approach to equilibrium would be governed by eigenvalue Xl+i- 
Here is the main result: 

Theorem 1. The spectral gap of the diffusion process (X(t))t>o * s given by r(/3, rj) = —9 where 9 
is the least negative solution to V(9;rj, f3) = with 

V(0;r?,/3) = -^D_ e {-(3)D'_ eh {JL) - D'_ e (-P)D_ e/v (fy (8) 

D u (z) the parabolic cylinder function with index v and argument z, and D' u (z) = ^D v {z). If f3 = 0, 
solving V = is equivalent to finding the roots of 

3/5 + I = (9) 

We also note that the result in (8) corresponds to taking the limit of the discrete queueing 
model with the scaling in (1) and then (7) examines what happens for large times. We show in 
Appendix D that (8) may also be obtained from the exact solution of the M/M/s + M queue. We 
have 

Proposition 2. The spectral gap in the discrete M/M/s + M model is the least negative solution 
to 

A(9) = F S (9)H S ^(9) - H S (9)F S ^(9) = 0, (10) 

where F n and H ri are the contour integrals in (199) and (201). For s — > 00, with po = 1 — P/y/s 
and 9 = 0(1), the roots of A(9) may be approximated by those of V in (8). This shows that the 
exchange of the limit in (1) and of large time is permissible in this particular case. 

Theorem 1 is an implicit description of the spectral gap, and it can be used to calculate r(/3,rj) 
numerically or asymptotically. For some values of /3 and 77 the spectral gap is shown in Figure 
1. We observe that the spectral gap decreases with /3 if rj > 1 and increases with /3 if rj < 1. 
If 7] = 1 the spectral gap is r(/3, 1) = 1 for all (3, since then the problem reduces to a standard 
Ornstein-Uhlenbeck process (see Section 3.2). This suggests that for systems with a large reneging 
rate, increasing the load (increasing po and decreasing f3) leads to shorter time scales for achieving 
equilibrium, while the opposite it true for small reneging rates. We also see that r increases as a 
function of 77. Later, we establish the monotonicity of r with /3, which is suggested by the numerical 
results in Figure 1 (see Section 7). 
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Figure 1: Spectral gap r(P,rj) for (3 G [—2,3] and n G (0,3]. 



To further substantiate our findings, we accompany the observations from Figure 1 by results 
for the spectral gap in various asymptotic regimes. In order to do so we assume that 77 — >■ (small 
abandonment rate). The asymptotics for 77 — > 00 can be obtained from the following important 
symmetry result, which we establish in Section 6. 

Proposition 3. For the density p there holds the symmetry relation 

p{x, t;x ;P,T}) = y/rj ■ p{-x^/rj, tr,; -x 0y /rj; ~P/y/rj, l/i]). (11) 

ForV in (8), 

V(0; v, P) = ~VV- V(9/ m l/ V ; -0/y/v) (12) 

and consequently, 

r(P,r j )= V -r(-P/ y /Tj,l/r ] ). (13) 

Next we give five different asymptotic results for r(P,rj) as rj — > 0. Then, from (13) we can 
immediately get results for 77—7-00. We shall consider five ranges of p, with P < 0, P ~ 0, 
< P < /?*, P w /3* and P > /3*. Here, /3* is the smallest positive solution to D'^ 2 , A {—P) = 0. We 
summarize in Table 1 the five cases and where the asymptotic result may be found. 

Proposition 4. For P < the spectral gap behaves asymptotically as r(P,rj) ~ rj with the correction 
term 



r(p, v )-r,~-2]§e 



ill 

/2tt" 

and hence r — n is exponentially small as i] — > 0. 



.si 

1 + Pe 2 / e 2 du 



00 



(14) 
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range of (3 


asymptotic result 


p < 
P > 0^ = 1.85722... 

o < p < A* 

P = 7VV = O(^rj) 
P- p* = tfl z W = 0(t/ 1/3 ) 


(14) 
(15)-(17) 

(18) 

(19) 
(23)-(25) 



Table 1: Five asymptotic regimes. 

Proposition 4 describes the part at the far right end of Figure 1, where r increases linearly with 
rj. Since P < 0, the diffusion process is mostly in the positive part of the state space, since the 
process has a positive drift for < x < —P/rj and an equilibrium point at x = —P/rj = \P\/rj S> 1. 
Hence, particularly when there is little reneging, one has to be far up in the state space before the 
process starts stabilizing. For the underlying queueing model, this scenario corresponds to large 
queues building up until enough customers renege so that the situation stabilizes. For this scenario, 
the spectral gap r = 0(rj) suggests large relaxation times. Note also from (6) that the steady-state 
distribution concentrates about x = —P/rj. Table 2 compares exact and asymptotic results for 
P = -l. 





P = 


-1 




V 


r (P, v) - 7 1 


(14) 




0.500 


2.50092-10" 2 


3.57325-10- 


-2 


0.400 


1.91877-10 -2 


2.48906-10' 


-2 


0.300 


1.16366-10 -2 


1.42105-10' 


-2 


0.200 


4.29814-10~ 3 


5.04257-10- 


-3 


0.100 


2.64792-10" 4 


2.92685-10- 


-4 


0.050 


1.32910-10 -6 


1.39448-10" 


-6 


0.025 


4.25017-10- 11 


4.47665-10" 
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Table 2: Results for p = -I. 

We next consider P positive and sufficiently large, where we obtain a very different result for r. 
Proposition 5. For P > p* = 1.85722 . . . and rj ->• 

r(P,rj)=r (P)+A(P)rj + O(r] 2 ), (15) 
where Tq{P) is defined implicitly as the minimal positive solution to 

D' ro (-P) = y / pyA-r D ro (-P). (16) 

The correction term is given by 



MP) - i P-VF^ D 

"2 p* - 4r ro 







-1 


. dp 


P=ro(0)- 





(17) 



where V(p,P) = D' p (- P) - D p (- P) y 7 /3 2 / 'A - p (so that V{r {P),P) = 0). We later show in Appendix 
E that A(P) > so that r(P,rj) — ro(/3) is positive for sufficiently small rj. 
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The equation (16) corresponds to the discrete part of the spectrum of the Halfin-Whitt diffusion 
with no reneging (i.e., with 77 = in (3)); see the discussion after Proposition 13. The various 
solution branches of (16) are demonstrated in Figure 2, where we plot the implicit function V(p, (3) = 
for j3,p > 0. 



P 



p = /3 2 /4 



Ps(/3) 
Pi 09) 



Figure 2: Solutions to vW 4 = pD p (-/3) = D' p (-f3). 

Proposition 5 applies to the flat part in Figure 1. Indeed, when (3 is large enough, the spectral 
gap is hardly influenced by 77. The diffusion process will spend most time below zero, near x = —f3. 
A likely queueing scenario would be that queues hardly ever build up, which makes the impact of 
reneging customers negligible. As the spectral gap is 0(1), we expect relaxation times that are 
0(1). While asymptotically r{(3,rj) ranges from t*o(/3) to 1, numerically this corresponds to the 
interval (/3 2 /4, 1) = (.86231, 1), which is quite small, leading to the flatness of the surface in Figure 
1 for /?>/?*. Table 3 compares exact and asymptotic results for (3 = 2. 

We next consider f3 > but with j3 < f3*, in which case (16) has no positive solutions. 

Proposition 6. For </?</?* = 1.85722 ... and 77 ->■ 0, 

KA V) = + V 2/3 \a \ (0 ^ + h (Kl + + °^ 4/3 )' ( 18 ) 

where ao = max{z : Ai(z) = 0} = —2.33810 . . . is the least negative root of the Airy function. 

From Proposition 6 we see that the asymptotic series now involves powers of 77 1 / 3 , which illus- 
trates the lack of analyticity of r(/3, 77) at 77 = 0. Now r(/3, 77) ~ /5 2 /4 and for 77 = the spectral gap 
is in fact exactly |/3 2 (see Section 3.1). Table 4 compares exact and asymptotic results for (3=1. 

When f3 becomes small, both (14) and (18) become invalid, as the correction terms become 
larger than the leading term. Then a separate analysis leads to the following result. 
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P 

r 


= 2 




t(B* 7l) 




0.5000 


0.98463 


0.93229 


0.2500 


0.97072 


0.93229 


0.1000 


0.95576 


0.93229 


0.0500 


0.94741 


0.93229 


0.0250 


0.94150 


0.93229 


0.0100 


0.93671 


0.93229 


0.0050 


0.93470 


0.93229 


0.0025 


0.93356 


0.93229 


0.0010 


0.93282 


0.93229 


±auie o. 


Results for /3 = 2. 




P 


= 1 


n 


r(j3,rj) 


(18) 


0.5000 


0.87510 


1.1778 


0.2500 


0.72686 


0.83452 


0.1000 


0.54242 


0.56732 


0.0500 


0.44074 


0.44990 


0.0250 


0.37193 


0.37593 


0.0100 


0.31673 


0.31836 


U.UUoU 


0.29217 


0.29306 


0.0025 


0.27664 


0.27713 


0.0010 


0.26450 


0.26472 


Table 4 


Results 


for 13 = 1. 



Proposition 7. Assume that (3 is small, such that (3 = 7y / ?7 = O(^rj) (f3 = •<-)■ 7 = 0) . Then 

r(p,rj)~riR(i), (19) 
where R is the minimal positive solution to 

jD R ( 7 ) = D 1+R (j). (20) 

Equation (20) has infinitely many positive solutions, whose existence follows from ODE theory, 
as discussed in (54)- (56). Note that using the relations (67) and (68) below, (20) is equivalent to 
Dr+i(j) + 2^(7) = or D' R {~f) + l^Djid) = 0. Also, R = is a solution for any 7. In Figure 3 
we illustrate the solution branches of (20) in the (7, R) plane, for R > 0. 

For R 7^ 0, (20) is also equivalent to Dr-i^) = 0, and then the solutions are precisely the 
positive eigenvalues of the reflected Ornstein-Uhlenbeck process (see the discussions in Section 3.2 
and surrounding (54)-(56) in Section 4). If 7 = we can use (66) to compute Dr_i(0) and its 
roots are R = 2, 4, 6, . . ., so that the spectral gap is r(/3, rj) ~ 2r/ if (3 = or (3 = o(y/rj). For certain 
special values of 7 7^ we can also get some of the eigenvalues more explicitly. For example, we 
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R 




-4 -3 -2-1 12 



Figure 3: A sketch of the solution branches of equation (20) for R > 0. 

know from (63) that D2 (z) is proportional to the Hermite polynomial He2(z) = z 2 — 1, so that 
Z?2(l) = 0, and then -R(l) = 3 is the minimal root of (20). We also have Z?2(— 1) = so that when 
7 = —1, R = 3 is a root of (20), but the minimal positive solution to (20) is R(— 1) ~ 1.3882 (which 
is illustrated in Figure 3). 

In Appendix A we establish: 

Proposition 8. For 7 — > ±00, R behaves as 

R-l ^Le~ 72/2 , 7^-00 (21) 

and 2 

R = \ + K|(t|) 2/3 [i + o(i)L 7^+oo. (22) 

Thus the approximation r ~ on the 7-scale in Proposition 7 asymptotically matches to the 
results in Propositions 4 and 6. In fact, in view of (22) the first two terms in (18) are a special 
case of (19), but this is not true for the third (0(rj)) term in (18). Note that r/7 2 /4 = /3 2 /4, 
jyy 2 / 3 = r? 2 / 3 /? 2 / 3 , and if (14) is expanded for small (3, r — r] agrees with r/(i? — 1) for 7 — > —00, in 
view of (21). 

The results in (21) and (22) are also consistent with Figure 1, which suggests that r(0, if) 
increases as a concave function of rj. The queueing counterpart is such that the load is one, and 
hence the reneging is necessary to alleviate the system. As rj becomes larger, more customers will 
leave the system, which reduces the queue lengths and, as seen from Proposition 7, shortens the 
relaxation times. 

Table 5 displays numerical results for 7 = 1,0, —1. 
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7 


= 1 


7 


= 


T 


-1 


1 


r(B, n) 

V" 5 1 J 


riRh) 


r(B, ri) 


riRh) 


r(B, w) 


riRh) 


0.50000 


0.81266 


1.50000 


0.65385 


1.00000 


0.54816 


0.69412 


0.25000 


0.53164 


0.75000 


0.38029 


0.50000 


0.29242 


0.34706 


0.10000 


0.24948 


0.30000 


0.16989 


0.20000 


0.12408 


0.13882 


0.05000 


0.13266 


0.15000 


0.08929 


0.10000 


0.06399 


0.06941 


0.02500 


0.06896 


0.07500 


0.04619 


0.05000 


0.03273 


0.03471 


0.01000 


0.02848 


0.03000 


0.01902 


0.02000 


0.01337 


0.01388 


0.00500 


0.01446 


0.01500 


0.00965 


0.01000 


0.00676 


0.00694 


0.00250 


0.00731 


0.00750 


0.00488 


0.00500 


0.00340 


0.00347 


0.00100 


0.00295 


0.00300 


0.00197 


0.00200 


0.00137 


0.00139 



Table 5: Results for 7 = 1,0, — 1. 



It remains to consider the case when 8 ~ 8*. Note that the correction term A(8) in (15) 
develops a singularity as 8 I 8*, since ro(/3) — > ro(/3*) = /3*/4. Also, the third (0(rj)) term in (18) 
becomes singular as 8 f P*, since by definition 8* satisfies D'^^i— /?*) = 0. Thus both (15) and 
(18) cease to be valid near 8 = 8* and we need a new expansion in this transition range. 

Proposition 9. For 8 ~ 8* such that 8 — 8*= r/ 1 ' 3 W with W = O(l), and rj — > 0, 



8* 



2/3 N 



with 



= \8*-^(^) 2l \(W) + 0{^ (23) 
where x is the maximal solution to 

Ai'(x)+(-|) 1/S -i-»'-Ai(x)=0 (24) 

L = i^hT) = 2 - 73875 • (25) 

IfW = {8 = 8*), then x(0) = max{z : Ai'(z) = 0} = -1.01870. . ., and as W ^ +co, \ -> +oo. 

As discussed in (57)-(59) in Section 4, the Sturm-Liouville ODE theory guarantees that there 
are infinitely many real solutions to (24). The solution branches of (24) are illustrated in Figure 4. 

We also note that if we order the roots of Ai(z) = as > oq > a\ > ■ ■ ■ and the roots of 
Ai'(z) = as > bo > b\ > ■ ■ ■ , these roots interlace as > &o > a o > b\ > a\ > • • • , and this 
fact can be used to establish more directly that (24) has infinitely many solution branches, for any 
fixed W. 

We have thus obtained the asymptotic connection between Propositions 5 and 6. Numerical 
results for the case 8 = 8* are given in Table 6. 

This concludes our asymptotic analysis of the spectral gap in Theorem 1. The proof of Theorem 
1 follows immediately from an explicit expression for the Laplace transform p of the transient 
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X 




Figure 4: A sketch of the solution branches of Ai'(x) + wAi(x) = 0. 







P = P* 




V 


r(/3,r/) 




(23) 


0.5000 


0.97803 


0.86231 


1.48841 


0.2500 


0.95673 


0.86231 


1.25673 


0.1000 


0.93129 


0.86231 


1.07644 


0.0500 


0.91493 


0.86231 


0.99721 


0.0250 


0.90139 


0.86231 


0.94729 


0.0100 


0.88770 


0.86231 


0.90845 


0.0050 


0.88016 


0.86231 


0.89138 


0.0025 


0.87462 


0.86231 


0.88062 


0.0010 


0.86966 


0.86231 


0.87225 



Table 6: Results for f3 = (3*. 

density, defined by 

/"OO 

p(x;6)= e- et p(x,t;x ;(3,r])dt, Re(0) > 0. (26) 
J o 

Defining the auxiliary function Ai by 

M(d- V ,(3) = ^D^ e (P)Di e/v (^) - DLMD.e/r,^), (27) 
we have the following result: 

Theorem 10. Consider xq < 0, with V in (8) and A4 in (27), and assume that Re(0) > 0. 
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(i) For x > 0, 



p( x - 9) = e^-^-^ D - 6i ~ X0 " P)D - e '^ ] (28) 

(ii) For x < xq, 

p(x;9) = e ^-^e\^-^ mD ~ e{ - X - ® 

V 27T 



(iii) For xq < x < 0, 



x ( D_ e (x + 0) + D_,(-x - ^(9^0) I ■ (29) 



V 27T 



x( flH ,(, + ffl + ^(-,-«^«). (30) 

The proof of Theorem 10 is presented in Section 6.2. Note that the results for xo > follow 
immediately from the symmetry relation (11). 

From (28)-(30) we see that singularities of p may arise either due to those of T(9) (which occur 
at 9 = 0, —1, —2, . . .) or from the zeros of V{9\ rj, /3) = 0. But in Appendix B we establish: 

Proposition 11. The singularities of (28)-(30) are precisely the solutions to V(0;r),(3) = 0. 

Hence, the large-time behavior of the diffusion process is dominated by the least negative zero 
of V, which gives the result on the spectral gap in Theorem 1. 

From (28)-(30), by evaluating the contour integral for the inversion of the Laplace transform 
p(x; 9), we can obtain a spectral expansion of the form 

oo 

p(x,t) =p(x,oo) + e^ {x °- x) e^ x o- r ' x2) Y, e ~ Xnt ^n(xo)^t(x), x > 0, (31) 



n=l 



and 



where 



p(x,t) =p{x,oo) + e^ x °- x ^e^- x ^^e~ x - t ^~{x )^-{x), x < 0, (32) 



and 



^(x) = VK rf - , x>0, (33) 

^n{x) = VK DX p { ~*J^\ x<0, (34) 

ju 1 n < mn ( & \ a* ^Hk3 ^ 

^ = A^ ( -^/H^)' = 89 e=-x n - (35) 
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The eigenfunctions ^ then satisfy the orthogonality relation 

r-O /-oo 



^n( x ) l l } m( x ) dx + {x)^m{ x )dx = 5(n,m). (36) 



o 

When r) = 1 we have A n = n, k n = D^) / \n\y/2Tr) , and D n (/3) = e -^/ 4 He n (/3) so that (36) 
reduces to 

/OO 1 ^ 
=He n (x + /3)He m (x + f3)e~^ x+ P )2 dx = 5(n, m). (37) 
-oo n!v2vr 

Note that the pole at = of (28)-(30) corresponds to the steady state behavior p{x,co), while 
the poles at 6 = —Xn and their residues lead to the decaying terms in (31) and (32). However, the 
spectral expansion does not yield any more insight than (28)-(30). 

3 Three special cases 

We shall now consider the three special cases of the diffusion process that arise by setting rj equal 
to zero, one and infinity. 

3.1 The Halfin-Whitt diffusion 

As rj —7- we end up with a process that behaves like a Brownian motion with drift above zero 
and like an Ornstein-Uhlenbeck process below zero. In [23] we have called this diffusion process 
the Halfin-Whitt diffusion, after Halfin and Whitt [17] who identified this process as a heavy-traffic 
limiting process for the GI/M/s system. The mean hitting time of the Halfin-Whitt diffusion was 
obtained in Maglaras and Zeevi [25]. Gamarnik and Goldberg [14] were the first to identify the 
spectral gap of the M/M/s system, asymptotically in the Halfin-Whitt regime. 

Theorem 12. (Gamarnik and Goldberg [14]) Let f3* = 1.85722... represent the smallest positive 
solution to D'^i^—p) = 0. The spectral gap of the M/M/s system, asymptotically in the Halfin- 
Whitt regime, is given by 

r ^' Uj " I r (/3), /?>/?*. (38) 

To establish Theorem 12, Gamarnik and Goldberg used the framework of Karlin and McGregor 
[20] for birth-death processes, and the result of Van Doom [9] on the spectral gap of the M/M/s 
system. In [14] the starting point is the discrete M/M/s model, and its spectral gap is then analyzed 
in the Halfin-Whitt regime (1). An alternative proof of Theorem 12 was given by the authors in 
[23] by deriving the expression for the Laplace transform p of the transient density in the diffusion 
limit, which shows that the limits of large time and (1) may be, in this case, interchanged. Below 
we summarize the main result in [23]. 

Proposition 13. For xq < and x > the Laplace transform of the transient density for the 
Halfin-Whitt diffusion with 7] = is 

- e^-f? 1 , ^ • (39) 

D-e(-P) VoTW/4 - D'_ e (-p)/D_ e {-P) 
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In Appendix C we show that (39) indeed follows by letting r\ — > + in (28). From (39) we see 
that there is a branch point at = —f3 2 /^, and this will lead to a continuous spectrum in the range 
lm(0) = and Re(0) < -/3 2 /4. There is a pole at = if (3 > 0, while if /3 < the pole is absent. 
Note that Do(—f3) = exp(— /3 2 /4). Other poles may occur at the roots of (16), which we studied 
analytically and numerically in [23] (see also Figure 2). 



3.2 Free-space OU process 

When rj = 1 it immediately follows from the process description that the diffusion process (X(t))t>o 
reduces to a free-space OU process, for which it is known that (with x > = max(x,xo), x < = 
min(x, xq)) 



p(x;6) 



e -M- x ^e^ { - x °-^T{e)D_ e {x > + P)D_ g (-x < - (3). 



(40) 



Indeed, this result also follows from Theorem 10 using the Wronskian identity in (69), which shows 
that 

V(0;r),P)^0-, r/^1. (41) 

Expression (40) is obtained for example in [34], in the context of the harmonic oscillator (see (46)- 
(50) below). Also, M -> (cf. (27)) as rj -> 1 and then (40) follows from (28)-(30). It is easy to 
invert the Laplace transform (40), as its poles are at zero and at all negative integers. Hence, (see, 
e.g., [21]) 



e \(xl-x 2 ) \p(x -x) °° r) t Xn + g)r> ( x + B) 
p(x, t; x ; /3, 1) = ±= ]T + P) ? n{ + P) e~ 

n=0 



2tt 



J?! 



(42) 



Here D n (z) = e z2//4 He n (z) where He n (-) is the nth Hermite polynomial. Alternatively, there is 
the closed-form expression 



p(x,t;x ;f3, 1) 



^/2tt V 1 - e' 



2/ 



exp 



(x + (3 - (x + /3)e-') 



t\2 



2(1 - e 



-2t> 



x e 



(43) 



3.3 Reflected OU process 

As rj —7- oo, the process will spend all its time below zero, and hence (X(t))t>o reduces to a reflected 
OU process (see Ward and Glynn [35], Linetsky [24] and Fricker et al. [12]). In this limit we have 
D-e/ v (P/Vv) -> A)(0) = 1 and y/vD'_ e/r) (P/^n) -> -/9/2, and then, using (67), 

V(0; V , (3) -> 9D^ (-P), r? ^ oc, (44) 
M{6; v, P) -> ^-i-e(^), " ^ oo, (45) 

which can be used to simplify (29) and (30). Then the solution agrees with that in Xie and Knessl 
[40, Eq. (4.3.2)]. 
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4 Schrodinger equations and spectral properties 



Here we give some basic background on spectral properties that are relevant to PDE's such as (2). 
In particular we show that the discreteness of the spectrum for any r\ > follows from classic 
results on the Schrodinger equation. 

We set p = e~ xt <p(x) where A is a spectral or eigenvalue parameter. Then (2) and (3) lead to 

(j)"{x) + {J3 + x)(j)'{x) + (A + l)(j>(x) = 0, x < 0, (46) 
4>"{x) + {13 + r]x)4>'{x) + (A + rj)<j)(x) = 0, x > 0, (47) 

and the interface conditions are </>(0~) = <^(0 + ) and </>'(0~) = <^'(0 + ). Furthermore, we can trans- 
form (46) and (47) into the self-adjoint form by setting 

\e-^e~^l^{x\ x<0, 
nX) -\e-W\-^ 2 /^(x), x>0, [ ' 

which leads to the Schrodinger equation 

- ip"(x) + V(x)tp(x) = Eip(x), -oo < x < oo, (49) 



where E and A are related by 



and the "potential" function V{x) is 



E=X + \, (50) 



V{x) = U(P + V xf + ^, x>0. (51) 

We also require the eigenfunctions ip(x) to satisfy ^(0 + ) = ip(0~), and ^'(0 + ) = i/j'(0~). Since the 
problem is defined over the entire real line, additional conditions must be imposed at x = ±oo, and 
most often it is required that ij)(x) € L 2 (— oo,oo), i.e., J*^^ \ip(x)\ 2 dx < oo. However, for parabolic 
and piecewise parabolic potentials, such as the one in (51), this condition is equivalent to simply 
rejecting solutions of (49) that have Gaussian growth as x — > +oo. 

We note that if the potential V(x) were exactly quadratic, say V(x) = x 2 /4, then the problem 
is just the quantum harmonic oscillator (or, for our application, the standard Ornstein-Uhlenbeck 
process), and then the eigenvalues are En = N +1/2 and the corresponding eigenfunctions are 
^n{x) = c^e~ x ^HeA^x), where cat is a normalizing constant and He]v(x) is the N th Hermite 
polynomial. Thus the spectrum is purely discrete for quadratic potentials. 

We can also view the differential equation in (49) as constituting a singular Sturm-Liouville 
boundary value problem. The study of such problems dates back to the work of Sturm in the 
nineteenth century, and they are discussed in detail in the books of Titchmarsh [34], Stakgold [29], 
Reid [28] and Coddington and Levinson [6, Chapters 7-12]. The problem in (49) is singular since 
it is defined over the infinite interval x £ (— oo, oo). 

Singular Sturm-Liouville problems are classified as either of limit circle or limit point type. For 
limit point problems the condition that the solution be square integrable is sufficient to determine 
it, while limit circle problems require a more explicit boundary condition at the singular point (s) 
(which are at x = ±oo for (49)). 
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Since Sturm-Liouville problems and Schrodinger equations are self-adjoint, their spectra are 
confined to the real axis. Singular problems may have both discrete and continuous spectra. How- 
ever, there is a general result, originally due to Weyl, with simplified proofs by Titchmarsh appearing 
in [32, 33] (see also the book [34]), that guarantees that (49) will have a purely discrete spectrum. 
This needs only the conditions that 

(i) V(x) be finite on finite intervals. 

(ii) V(x) — > +oo as x — > ±oo. 

Our potential in (51) clearly satisfies these conditions and thus has a purely discrete spectrum, for 
any rj > 0. The fact that V(x) has a jump discontinuity at x = does not affect the spectrum; 
it only means that some jump conditions must be specified at x = 0. However, if rj = 0, then the 
potential does not grow at x = +oo, and then in fact, as we discussed in [23], the problem has a 
continuous spectrum in the range A > /3 2 /4 (E > /3 2 /4 + 1/2), and may also have any number of 
discrete eigenvalues, depending on the value of f3. Much of the asymptotic work here assumes that 
rj — > + , so we are looking at a very singular limit where the discrete spectrum begins to resemble 
a continuous one, in certain ranges of A. 

For the problem in (51) the smallest eigenvalue is Eq = 1/2 (thus Ao = 0) with the corresponding 
eigenfunction being the piecewise Gaussian 



V>oOr) 



e -/3x/2 e -x 2 /4 ! X<Q 
e -/3x/2 e -^ 2 /4 ) x > o 



and this corresponds to the steady state distribution in our model. 
Given the discrete spectrum we order the eigenvalues En as 



Eq < E\ < Z?2 < • • • < E^ < . . . 



(52) 



with A^v = En — 1/2. By general results for Sturm-Liouville problems the sequence {En} satisfies 
En — > oo as N — > oo. Also, for every eigenvalue there is only one linearly independent eigenfunction, 
so all eigenvalues are simple. This can be shown directly from (49), for if tp(x) and ip(x) corre- 
sponded to the same eigenvalue E, then ^(x)^" (x) — ^"(x)tl>(x) = ^(x)i/j' (x) — ^ / (x) , 0(x) = 

so that the Wronskian ipip' — ip'tp is a constant, say c. But by letting x —> ±oo we conclude that 
c = so that 



d_ 

dx 



ip(x) 
tp(x) 







(53) 



and thus ip must be a multiple of ip. Since V(x) in (51) js discontinuous at x = 0, we can use (53| 
to conclude that i(j(x) = ktp(x) for x > and ip{x) = k'ip(x) for x < 0, but then continuity of 
at x = forces k = k! . 

There are two other singular Sturm-Liouville problems that are relevant to the analysis here. 
First consider 



- V"(X) + 
with the boundary condition 



*(x) = m(x), o < x < oo 



*'(0) + ^tf(O) = o. 



(54) 
(55) 
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This problem has a regular point at X = 0, where a standard boundary condition is applied, and 
a singular point at X = oo, where we require that ^(X) 6 L 2 (0,oo). This is a singular problem of 
limit point type at X = oo which may be explicitly solved in terms of parabolic cylinder functions, 
with ^(X) = Dg(X + 7). Then (55) leads to the eigenvalue condition 

D'~( 7 ) + ^(7) = 0. (56) 

The results of Weyl and Titchmarsh again guarantee that the problem has a purely discrete 
spectrum and thus an infinite sequence of eigenvalues {E^}. Also, E = is the lowest eigenvalue 
with ^o(X) = e~^- x+ ^ / 4 . Note that (56) is essentially the same as equation (20) in Proposition 7, 
since D'^il) = \l^E^ ~ ^E+iil)- Thus the existence of infinitely many real solutions to (20) 
follows from Sturm-Liouville ODE theory, though in the next section we shall also show that it 
follows from the oscillatory nature of the parabolic cylinder functions, as functions of their index E. 

Another singular Sturm-Liouville problem is 

— i/j"(x) + xip(x) = Etp(x), < x < 00 (57) 

with the boundary conditions 

<//(0) +w0(O) = (58) 

and ijj(x) S L 2 (0, 00). Here a; is a real parameter. Again, since V(x) = x grows linearly as x — > 00 
and x = is a regular point, we have a purely discrete spectrum. But (57) is related to the Airy 
equation, with solutions proportional to Ai(x — E), and the eigenvalues are determined by (58), 
hence 

Ai'(-E)+ujAi(-E) =0. (59) 

But (59) is equivalent to equation (24) in Proposition 9, so again ODE theory can be used to 
establish the existence of infinitely many solutions. Note also that if 00 = the eigenvalues are the 
roots of Ai'(-), while if oo = 00 the eigenvalues are the roots of Ai(-). 

To summarize we have given some basic background on Sturm-Liouville theory, and on Schrodinger 
equations and their spectral properties, that are useful in the present investigations. In particular 
this theory guarantees infinitely many discrete solutions to the equations that arise in Propositions 6 
and 7, and in Theorem 1. 



5 Parabolic cylinder functions and their properties 

The parabolic cylinder equation is the second order ordinary differential equation 

y"(z)+(p+^-^)y(z) = 0, (60) 

where z is a complex variable and p is a parameter. Since (60) has no singular points (except at 
z = 00) its solutions are entire functions of z (see [6]). 

One solution of (60) is denoted by D p (z), which is called a parabolic cylinder function of order p, 
and it is defined by the integral representation 

D p {z) = -^=e z2/i [ u p e- zu e u2/2 du. (61) 

W2lT J Br 
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Here Br is a vertical Bromwich contour on which Re(it) > 0, and the branch of u p will be defined 
by u p = |u| p e iarg (") where —it < arg(tt) < it. Then the integrand in the right-hand side of (61) is 
analytic exterior to the branch cut where Im(u) = and Re (it) < 0. The function D p (—z) provides 
a second linearly independent solution to (60), so that 

Cl D p (z) + c 2 D p (-z) (62) 

is the general solution, with c\ and C2 being complex constants. When p = 0, 1, 2, . . . is a non- 
negative integer we can obtain D p (z) in a closed form as 

D p {z) = D N (z) = e- z2/A Ue N (z); N = 0, 1, 2, . . . (63) 

where Hejv( - ) is the ./V th Hermite polynomial. Note that if p = N the integrand in (61) becomes an 
entire function of u. Here we use the notation He(-) for the Hermite polynomials, so that Heo(z) = 1, 
Hei(z) = z and in general HeAr(z) ~ z N as z — > oo. When p = N, D^{—z) = (—l) N D^(z) and 
another linearly independent solution must be used in (62), but we shall not need it in the present 
analysis. The function D p (z) is real valued when z andp are real. As discussed in [31], D p (z) is an 
entire function of both p and z, and indeed one can easily compute derivatives of all orders from 
the integral representation in (61). For example, we have 

BtW = = / (f - u) u>e-°«e^ d u (64) 



and 

9 -DJz) = -^=e z2/i I u p {\ogu)e- zu e u2/2 du. (65) 



dp py ' 



Br 



In (65) logtt is real for u real and positive, and analytic exterior to the cut Im(n) = 0, Re(n) < 0. 
If z = the integrals in (61) and (64) may be expressed in terms of the Gamma function, with 

^ (0) = f7T^y d p {0) = r (^) " ( } 



Since T(z) has simple poles at z = 0,-1,-2,... it follows that D p (0) has simple zeros at p = 
1,3,5,..., while D' p (0) has simple zeros at p = 0, 2, 4, . . . . This also shows that the functions in 
(66), as functions of p, tend to oscillate for p > 0, but have one sign for p < (actually, for all 
p < 1 for D p (0)). Also, in view of the growth of T(z) as z — > +oo, the functions in (66) decay very 
rapidly for p — > — oo. 

Using (61) and (64) we can easily derive the recurrence relations 

D' p {z) + \zD p (z) - pDp_i(z) = 0, (67) 

D' p {z) - \zD p (z) + D p+l {z) = 0, (68) 

which we shall use in the present analysis. 

The Wronskian of D p (z) and D p (—z) is defined as D p (z)D' p (—z) + D' p {z)D p {—z) and it has a 
very simple form, with 

D p (z)D' p (-z) + D' p {z)D p {-z) = -^- y 
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which vanishes if p = 0,1, 2,... 
neously. For if D p (z*) = D' p (z* 



We also note that D p {z) and D ' Jz) cannot both vanish simulta- 
= for some p and z*, then (60) shows that D''{z*) = 0. Then 



repeated differentiation of (60) would show that all derivatives of D p (z) vanish at z 



Then we 



could expand D p (z) in Taylor series about z = z* to conclude that D p {z) = in some neighborhood 
of z = z*. But since D p (z) is an entire function of z this would imply that D p (z) is identically zero, 
which is clearly not the case. 

To better understand the behavior of these special functions, many asymptotic formulas have 
been derived for z and/or p large. We summarize some of these below, since they are used to 
establish our main results. First, for z large and positive, we have (see for example [16, p. 1093]) 



D p (z) 



7i 



1 



Pip - 1) 

2z 2 



+ 0{z 



oo, I arg(z)| < 37r/4. 



(70) 



A more general result, which allows z to be negative, is the following (see for example [16, p. 1094]): 



D p (z) 



1 - P[P - 1] + O(z^) 



2z 2 



(71) 



j^^/V- 1 [1 + 0(z" 2 )] , z -> oo, \ < arg(.) < ^. 



X-p) 

Here we let z = \z\e % ar §( z ) where | • | denotes the complex modulus. 

The limit of z large and negative corresponds to setting arg(z) = it in (71), and then the leading 
term becomes 



D p {z) 



\-z) 



-p-i 



-oo, 



T(-p) 

which holds as long as we are away from the zeros of l/T(—p), which occur at p = 0, 1, 2, 



(72) 



If 



N= 0, 1,2, . . . we have D p (z) = e z / 4 Hejv(z) and then (70) holds for all values of arg(z) 



P = 

In some of the analysis that follows we will need to consider cases where p is not exactly a 
non-negative integer, but is very close to one. For p ^ iVthe first series in (71) has Gaussian decay 
as z — > — oo (0(e~ z / 4 )), while the second series in (71) (cf. also (72)) has Gaussian growth as 
z — > — oo. But if p is very close to an integer these two terms may be of comparable magnitude. 
For example if p = e is small T(—p) ~ — 1/e and if e — > and z —> — oo in such a way that e z l 2 e 
is 0(1), then the two parts of (71) are roughly comparable. 

The asymptotic results in (70)— (72) follow easily by expanding (61), using techniques for the 
asymptotic evaluation of integrals, such as the saddle point method and singularity analysis. Gen- 
eral references for such methods are the books of Bleistein and Handelsman [5], Wong [39], Sz- 
pankowski [30], and Flajolet and Sedgewick [11]. The integrand in (61) has a saddle point at 
u = z and a branch point at u = 0, and one of these (or both) determines the asymptotic behavior 
of D p (z) as z — > oo, for any direction arg(z) in the complex plane. The saddle leads to (70) and 
the first part of (71), while the branch point leads to (72) and the second part of (71). 

We next consider a fixed (real) z and expand D p (z) in the limits of p — > dboo. Then 



D P (z) ~ -j= exp 



P -i / \ P . 

- lOg(-p) - - + Zy 



P 



D p (z) = V2 



exp 



log p 



cos (ptt - Vp z ) + 0(p x ) 



-oo, 



p 



+oo. 



(73) 



(74) 
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From (74) we see faster than exponential growth with p, coupled with oscillations, in view of the 
trigonometric factor. Near a zero of the cosine the 0(p _1 ) error term may become important, and 
it may be also explicitly obtained. Thus D p (z) has an infinite number of zeros as p increases toward 
+00, for any fixed real z. This is in sharp contrast to viewing D p (z) for a fixed p as a function of z, 
in which case it has at most finitely many zeros. It is known [1, p. 696] that D p (z) has no zeros in 
the range p + 1/2 < z 2 /4. Also, (74) shows that the large zeros can be estimated by 



p = 2M+ 1 + 2-V2M+ 1 + O(l) 

7T 



M- 



00. 



(75) 



The results in (73) and (74) may be obtained, for example, by expanding the integral in (61). When 
p — > —00 the asymptotics are governed by a single saddle point at u = y/—p, while as p — > +00 two 
saddle points, at u = ±i*Jp, contribute. 

We next consider asymptotic limits where p and z are simultaneously large, restricting ourselves 
to real p and z. The results quoted below are taken out of Abramowitz and Stegun [1], where 
we note that in [1] the results are given for the function U(a,z), which is related to D p (z) by 
U(a, z) = D_ a _i/ 2 {z), or D p (z) = U(—p — 1/2, z). A complete summary of the asymptotics of U is 
also given in Temme [31]. 

When p + 1/2 > and z 2 — 4p is large and positive, the so-called Darwin's expansions apply, 
where [1, pp. 689-690] 



D p (z) 



Vr(p + 1) 
(2tt)V4 



x exp 



4p-2) 



-1/4 



1 + o[ (z 2 -4 P y 3/2 



1 



4p - 2 + [p + - lo 



4p-2 



(76) 



and for p + 1/2 < and z 2 — 4p large and positive 



D p (z) 



(2k) 



1/4 



x exp 



4p-2) 



-1/4 



Ap 



1+ 0[(z 2 -4p) 



-3/2 



(77) 



The results in (76) and (77) apply for \z\ — > 00 and p — > ±00 with z 2 /(4|p|) fixed, and are uniform 
in the interval z 2 /(4|p|) G [1 + e, 00) for any e > 0. The expressions in (76) and (77) are more 
uniform than (70), (72) and (73), and contain these as special cases. 

Now consider z and p large with z 2 — 4p — > —00. For p > the appropriate expansion is now [1, 
p. 690] 



2^r( P + i) 2 1/4 

p(z) ~ (2tt)V4 ( 4 P + 2 ~ z ) 
7rp z 



(78) 



X cos 



p + - J sin 



-1 
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Here sin _1 (-) G (— tt/2, it/2). The result in (78) is uniform for p — > oo and z — > ±00 for intervals 
where z 2 /{Ap) G [0, 1 — e] for any e > 0, except if we are at or near a zero of the cosine function, 
in which case the correction term(s) to (78) must be considered. 

For a fixed large z, p has to increase past approximately z 2 /A in order for the zeros of D p (z) to 
become evident. The expansions in (76) and (78) develop non-uniformities when z 2 /{Ap) ~ 1 and 
there yet other expansions apply. The following result [1, p. 689] is more uniform and applies for 
all z 2 /(4p) G [0,oo], as long as p and \z\ are large: 



DM ~ 2^T (E±i) (4 P + 2) V (^J-) Ai ((4p + 2f 3 r) 



and 



V47+2' 



2/3 



(79) 
(80) 



geVe^i-gCosh-^oj ,t>i. 



(82) 



Here Ai(-) is the Airy function, which has the following asymptotic behaviors as z — > ±00 (see [1, 
p. 448]) 

Ai(z) ^z" 1 / 4 exp ( -lz 3 A , z -> +00 (83) 
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Ai(z) ~ ^(-^)" 1/4 (^("^) 3/2 + £ ) , * "> "oo- (84) 



For p — > 00 with a fixed r > 0, we can simplify (79) by using (83) to approximate the Airy function, 
and then we obtain (76) as a special case, when £ > 1. 

Similarly, for £ < 1, we can use (84) to approximate the Airy function in (79) and then (79) 
reduces to (78), up to a Stirling approximation for the Gamma factors. When £ = 1 (r = 0) or 
£ 1 we can also simplify (79) considerably. Suppose that p — > 00 and r — > in such a way that 
p 2 I^T is fixed. Setting 

(85, 

with 5 fixed as z — > 00, we have (4p + 2) 2 / 3 r ~ 5 and then (79) simplifies to 

D p (z)-e-^Q P V^Q 1/3 M(S), (86) 

We shall use (86) to establish Propositions 6 and 9. 

We have thus summarized the various "uniform" asymptotic approximations to D p (z), where 
both z and p become large. Despite the seeming complexity of these results, they are easily obtained 
from (61) via the saddle point method. Indeed, setting u = zv (with z > 0) in (61) leads to 

D p (z) = Z * \ e^'^dv, $(v) = l-\ogv+ V --v. (87) 
W2ir Jbx z 2 
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The saddle point equation is <&'(v) = which is the quadratic equation 



v 2 -v + ^ = v = v± = -\l± v/l - 4pz- 2 l . (88) 

To expand the integral in (87) for z — > oo, p — >• oo with p/z 2 fixed we find that for 4p/z 2 < 1 the 
real saddle at v = v+ determines the asymptotic behavior, and we ultimately obtain (76), up to a 
Stirling approximation of T(p + 1) and the equivalence z 2 — 4p — 2 ~ z 2 — 4p. In contrast, when 
4p/z 2 > 1 two complex saddles, at i\J 4pz~ 2 — l] , contribute to the asymptotics and then 

we obtain (78). The transition range in (86) corresponds to 4pz~ 2 ~ 1 and then the two saddles 
coalesce to form a higher order saddle, and such transitions invariably involve Airy functions (see 
Chapter 9 in [5]). 

We have discussed here only approximations to D p (z), but some of our main results involve 
also the derivative D' p (z) (see, for example, Theorem 1). Its asymptotics follow from the integral 
in (64), but the same results can be obtained by formally differentiating the results for D p (z), as 
in this case term by term differentiation of the asymptotic series is permissible. For example, the 
logarithmic derivatives of D p (z) and Ai(z) satisfy 



D ' P ( Z ) z , P , P(P~ 1) , n( -5\ 
Djz)=-2 + z + —^ + °( Z 

Ai'(z) r- 1 



(89) 



yfz - — + [ z~ 5/2 ) , z -)■ +oo 



Ai(z) Az 
and later we shall make use of these results. 

6 Proofs 

6.1 Proof of Proposition 3 

Here we establish the symmetry relations (11)-(13). These may be obtained without solving ex- 
plicitly for p(x, t). Consider the problem in (2) with p = p(x, t; xq; (3, rj) and set 

x = ~, t = -, P = -P'^Tj, x = -^L (90) 



with 

Then (2) becomes 



p(x, t; x ; /3, V ) = R(x',t'; x' ; 0, r/). (91) 
■qR t > = 7]R x r x > + T]j3 R x i + x R x i + R, x > 0, 

r)R t , = r]R x , x , + r)$ R x > + ri[x'R x , + R], x' < 0, (92) 
where R x * = dR/dx' and the initial condition becomes 

R\t'=o = 5 (r^7=f-) = ^ x 'o - x ') = - x 'o)' ( 93 ) 
where we have used the scaling law of the delta function. 
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Dividing (92) by rj and setting R = y/rjR we see that p and R satisfy the same problem in the 
x',t' variables, with 77 replaced by 1/rj. Hence 

-^-R(x', t'; x' ;/3', rj[) = p(x', t'; x' ; 0, 1/r/) (94) 
V 7 ? 

and then (91), with (94) and (90), is equivalent to (11). 

To establish (12) we again replace /3 by —f3'y/rj and also set 9 = r/9' . Then using (8) we obtain 



d 



1 d 



V(9; V ,/3) = -^D_ v9l (/3'^)[-—D_ ei (-p')\+——[D_ ve ,(P'^ 

= -y/rjV{#\ l/ v , p) = -y/TjV(9/ri; 1/r?, -0/Jrj), (95) 

which establishes (12). Then since by definition r = r(/3,rj) is the minimal negative root of V = 
0, the right-hand side of (95) has a root where —9/t] = r(— /3/y/fj, I/77) and then (13) follows 
immediately. 

6.2 Proof of Theorem 10 

We let p(x; 9) = J* °° e~ et p(x,t)dt and note that p will be analytic in the right half-plane Re(0) > 0. 
If p satisfies (2) its Laplace transform satisfies 

9p(x; 9) - 5(x - s„) = ~[a(x)p(x; 9)] + , (96) 

where 

- ^ [a(x)p(x, 9)\ -j {x + ^ J[x . Q) + p{x . Ql x<Q (97) 

Assume that xq < and x > 0, so that 5{x — xq) = 0. By writing p = e~ r,x ' 2 / 4: e~l 3x ^ 2 v(x; 9), (96) 
reduces to the differential equation 



11 

v 



{x-9) + {\r,-9-\{xi 1 + l3) 2 )v{x-9) = 0, x > 0, (98) 



where v' = dv/dx. This is the parabolic cylinder equation (Erdelyi [10], p. 116) and as we discussed 
in Section 5 two linearly independent solutions (at least for Re(#) > 0) are given by 

D-0 /v C-^), D_ eh {^). (99) 

Thus we write 

p(x; 9) = l4 (9)e-^/*e-W 2 D_ e/r! (^), x > 0, (100) 

as the second solution in (99) must be rejected due to its Gaussian growth (see (72)). For x < we 
must use the second expression in (97), and then solve (98) with 77 = 1. Now we reject solutions 
with Gaussian growth as x — > —00, so that for x < xq the appropriate solution to the second 
equation in (97) (using (98) with r] = 1) is 

§{x-9) = j 1 (9)e-* x2 e-^ x D„ e (-/3-x), x < x < 0. (101) 
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But in the range xq < x < the solution will involve both of the parabolic cylinder functions 

D_q(—/3 — x) and D_g((3 + x), hence 

p(x; 9) = e -i*V> [ l2 {9)D^ e {-p - x) + iz{9)D-e{P + x)] . (102) 

The functions 7j(0) are determined from continuity conditions at x = and x = xq (cf. (4) and 
(5)). Continuity of p and -jj^p at x = leads to 

j 2 D-g(-P) + -ftD-oiP) = 7iD_ e/v (^), (103) 
- l2 DL e (-(3)+ l3 D'_ e ((3) = lA ^D'_ eh {JL). (104) 

Continuity of p at x = xq yields 

TiD-e(-P ~ x o) =i2D- e {-P - x ) + j 3 D_ e (p + x ), (105) 

and the jump condition of -^p at x = Xq, i.e., 

p{x+;6)-p{xQ-9) = - f ° 6(x-x )dx = -1, (106) 

leads to 

-e^es^o = - l2 D'_ e (-P - x ) + j s D'_ e (P + x ) + iiD'_ e {-P - x ). (107) 

Here we used (101) to compute p(xq), (102) to compute p(xq), multiplied (106) by e^ x °e^ x °, 
and used also (105). Equations (103)-(105) and (107) give a 4 x 4 linear system for the jj, whose 
solution leads to Theorem 10. 

We note that the calculations assumed that Re(#) > 0. If Re(#) < the 4x4 system may 
become singular, and in fact this occurs when 9 = and at the eigenvalues —9 = Atv, N > 1. 
Theorem 10 thus gives the Laplace transform p for Re(#) > 0, and then the expression can be 
analytically continued to the left half-plane, since we know how to continue the parabolic cylinder 
functions, which are entire functions of 9. After the continuation, locating the singularities in 
(28)-(30) in the range Re(9) < can be used, for example, to obtain the spectral representation in 
(31). 

6.3 General considerations for establishing Propositions 4-9 

Here we discuss some general principles about solving V = in Theorem 1, for the minimal root 
r(/3,r]), in the limit of rj — > + . As discussed in Section 4, general results for Schrodinger equations 
and Sturm-Liouville problems show that the roots of V = are all on the real axis, and that the 
sequence of roots (or eigenvalues) —9 n = A at satisfies A at — > oo as N — > oo, for any fixed j3 and 
r] > 0. Also, we know from Section 4 that the roots are all simple, and thus dV/d9 ^ when 
9 = —Xn- 

Consider r as a function of /3 and r]. Then V(— r(/3, rj); T], f3) =0 and by implicit differentiation 
we obtain 

dV dr dV , s 

89 e=-r dp dp 
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and since dV /d9\e=- r / we can use this relation to compute dr/dfi (for any /3 and any 77 > 0). 
By taking higher order derivatives of V = with respect to /3, a similar argument shows that r has 
derivatives of all orders with respect to j3. Also, by differentiating V = implicitly with respect to 
77, we conclude that r has derivatives of all orders with respect to 77, for any 77 > 0. Thus r(/3, 77) is 
infinitely smooth for all real /3 and for 77 > 0, and this is true for the higher roots also. Note that 
since D p (z) is an entire function of both p and z (see [31]), V is an entire function of 9 and (3, and 
real analytic for 77 > 0. However, the limit 77 — > + is quite singular, as we shall show. 

The above discussion shows that the roots of V = vary smoothly with /3 and 77, and a root 
cannot simply appear/disappear, say at some critical value 77 c . Thus for 77 — > + the roots have to 
lie in some range(s) of 9. In Section 5 we gave detailed results of the different asymptotic expansions 
of the parabolic cylinder functions D p (z), for different ranges of p, z. Applying these results to the 
equation V = 0, the function V can be approximated by simpler functions in the limit of 77 — > + , 
but these approximations are different in different ranges. 

Different expansions arise for the ranges —6 = 0(77), < —9 < /3 2 /4, —9 « /3 2 /4 and —9 > /3 2 /4. 
Depending on the value of /3, we shall need to consider different ranges of —9 in order to locate the 
minimal root. 

In what follows we shall use the following principle: suppose the equation F(u, e) = has roots 
Uj = Uj(e) which depend on the small parameter e, and these roots are smooth functions of e. Also, 
suppose that F(u, e) is an analytic function of both u and e, with an expansion of the form 

F(u,e) = F (u) + eF 1 (u) + 0(e 2 ). (109) 

Then if Fq(u) has a simple root at then F(u,e) has a root close to u* for e — > 0. The same 
conclusion holds if F is not analytic in e, but has an asymptotic expansion of the form (109), where 
the expansion holds uniformly on some (finite) u interval that contains tt*. In our case —9 plays 
the role of u and e will correspond to 77 or a fractional power of 77, such as ^Jrj or 77 1 / 3 . 

6.4 Proof of Proposition 7 

For (3 = 0(y/rf) and 77 — > + , we shall show that V = has solutions in the range —9 = 0(rj). We 
use the facts that Dq(0) = 1 and D' (x) ~ —x/2 as x — > 0. We let 77 — >• 00 and from (8) we obtain 

V(9;ri,l3)^^D-e(-P)-D'_e(-P), 77^00. (110) 

The error term in (110) is uniformly 0(?] _1 ) on finite 9 intervals, and we note that the right-hand 
side of (110) is an entire function of 9, as will be the error terms. 

Then the symmetry relation for V in Proposition 3 implies that if we scale 9 = r]S and /3 = ^y^/rj 
we obtain 

V(0;77, / 3)~-07[ J D_ 5 (7)| + J D-s(7)] ) V -> 0, (111) 

with an error term that is uniformly 0(r]) on finite intervals of S and 7. Proposition 7 follows upon 
setting S = —R and using the identity Di-s(j) + D'_ s (j) = ^~fD-s(j) (see [16], p. 1066). 

Finally, we show that for f3 = 0(y/rf) and 77 — > + , there can be no roots of V = in (8), other 
than 9 = 0, in the range 9 = 0(77). Consider the scales 9 = Q(rj M ) for M > 1, setting 9 = 9^r\ M . 
Then we have 

V = -^D_ e ^ M (-y/rjj)D _ m m_i(7) " D' M {-^)D_ e ^M-iin). (112) 
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Then both 9*rj M and 6»* i] M 1 are small and expanding the parabolic cylinder functions in (112) in 
Taylor series in both index and argument yields 



V 



D' 



.(7) 



D' {- y /Tj 1 ) + O(9^ M ) £L m m( 7 ) 



M M - 1/2 ^K(7) + i^(7)] | p= Jl + o(l)] 
M M ~ 1/2 ^-i(7)[l + o(l)], 



(113) 



where the error term is o(l) for 77 — > + and this is uniform in finite 9* intervals. But then we 
conclude that 6* = 0, contradicting our assumption that there is a root in the range 9 = Q(rj M ) 
for M > 1. To obtain the last expression in (113) we also used the recurrence (67). 

The calculation that led to (113) only used the fact that 9 and 6/rj are both small. Indeed, 
for any 9 = 0(77) we obtain (113) with 9*rj M replaced by 9. Thus if there is a root in any range 
where 9 = o(rj), we again conclude that 9 = 0. Hence, there can be no roots in ranges where 



etc. 



•l°g(l/»?) ; ' V log log(l/»7) 

6.5 Proof of Proposition 5 

Here we establish Proposition 5. We represent the parabolic cylinder function (see also (61) and 
(87)) by the contour integral 



-o/v 



i\/2ir Jc 



t -6h e -^ e tV2 dt 



e P 2 /4v 



W27T 



1 
Tj 2 ') 2 



cxp 



dz, 



(114) 



where (j) = ~ fiz — 9logz and C is a vertical Bromwich contour with Re(t) (or Re(z)) larger 
than zero. Now assume that 9 + j3 2 /4 > and 77 — > + . From the discussion in Section 5 we can use 
(77) (or the more uniform result in (79)) to approximate D_ e / ri (f3/^/rj) in this range. We briefly 
derive the asymptotic formula below. 

The integrand in (114) has a saddle point where <fi'(z) = 0, which occurs at 



+ + 46> 



lis) 



For /3 2 + 49 > this saddle point lies on the real axis and the directions of steepest descent are 
arg(z — z*) = ±7r/2. Then expanding <p in Taylor series about z = z* and noting that 



e H**mh exp 



1 

2rj 



<t>"(z*;P,9)(z-z*f 



dz = exp 








7] \ 



fiz* - 6>logz* 



(116) 



we obtain 



D_ 



oh 



r) 2 v 1 + — 



-1/2 



exp 



\{\zl - Pz* ~ & log 



(117) 
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In view of (64), the integral representation of D'_Q,((3/y/rj) corresponds to multiplying the inte- 
grand in (114) by {13/2 — z)j ^Jr\ (after scaling t = z/^/rj). Then again applying the saddle point 
method we see that the leading term in the expansion of ^/rjD'_ s ^(f3 / \/rf) is essentially the same 
as (117), with the additional factor /3/2 — z*, and thus 



vW_ e/Tj (/V^) 

D-e/niP/Vv) 



P 



(118) 



In (118) we divided by D_Q/ v ((3/y/rj), since this is nonzero in the range 9 + /3 2 /4 > 0, as discussed 
in Section 5. Using (118) we see that V = leads to (after dividing by D_ e / r] {j3 / \frf)) 



Z>_ e (-/3)vW4 + - D'_ (-P) + 0{ V ) = 0. 



(119) 



For ))->0we obtain the limiting equation in (16). Up to now the calculation did not distinguish 
between or ft > /3*, but in the former case (16) (or (119) with rj = 0) has no roots (other 

than 6 = if j3 > 0). By computing the first correction term in (117), which is 0{r}) relative to the 
leading term, we can refine (118) to 



(120) 



Here the error is uniform for 77 — ?- with \J /3 2 /4 + 9 > e > 0. However, the asymptotics break 
down as /3 2 /4 + # — > 0, and a separate analysis is needed for locating the roots of V = in the 
range — 9 « /3 2 /4, which we carry out in the proofs of Propositions 6 and 9. With (120) we can 
improve (119) to 



^( P --^T9) +0 ^) 



P 2 + 49 V2 



0. 



(121) 



Then the leading term in (15) follows by letting r\ 
dividing (121) by r/ and then letting 77 — ^ and — 1 

1 



in (121), and the correction term follows by 
-> tq simultaneously, noting also that 



lim 



D-e{-P)y/p/4 + - D>_ e (-P) 



lim 

r/— >0,— 6— >rro 



-v(-M) 
v 



dV 



■ lim 

=r r?-»0 



r(P,ri)-ro(P) 
V 



(122) 



where by definition in (15), A(/3) = lim r? _ s> o[ r (/3, f]) — ro(/3)]/T]. 



6.6 Proof of Proposition 4 

For {3 < we analyze the range 9 = 0(r/), and we shall see that there are roots in this range. Note 
that for /3 > the calculation in the previous subsection (since (118) applies for 9 = 0(rj)) shows 
that the only root in this range is at 9 = 0. Since we are examining ranges where 9 is small, we 
again use the fact that D_$(—f3) is an entire function of 9, and hence by Taylor series, noting that 
D (-P) = e"^ 2 / 4 and D' {~P) = ^e^ 2 / 4 , we have 



1 jjD-Bi-P) l_ 



l-9TZ(P) + 0(9 2 ) 



(123) 
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where (with z = —/3) 



d D' p (z) 



dp D p (z) 



e 2 2 /2 r i 
e ' / 1 t 2 /2-zt 



p=0 



iy/2ir Jc t 



=/3 2 /2 



dt 



- u2 l 2 du. 



(124) 



We then rewrite V = as 

i w D -»(-P) 



1 



l-9K{[3) + 0{9 2 ) 



D_ e/v (f3/^j), 



and rearranging terms we obtain 



d 



^=n{P)D-e h {^) + 0{9\ 



(125) 



(126) 



Note that 9 = is a solution to (126). Now consider 9 < with —0 = — t/r/ = 0(rj). For /3 < and 
?7 — > 0, (72) shows that 



D q (P/rf)=D q (-\P\/y/fj) 



r( 



(127) 



as long as<77^0,l,2,.... Using a similar formula for D' q ((3/^/rj) we then have 

^(/Vv^?) /8 



9^0,1,2, 



(128) 



But then (125) cannot be satisfied asymptotically. We conclude that if V = has roots in the range 
9 = 0{rf) they must occur where —0/rj = q fa 1, 2, . . . (in addition to the root at 9 = 0). To obtain 
the minimal root we examine the range where q ~ 1, thus setting q = 1 + e where e = e(rj) — > 
(which corresponds to —9 = rj + r/e(r/)). Then (126) is equivalent to the asymptotic relation 



Vv 



d_ 

dp 



e 4r > D 



l+e\ 



^lZ(P)D 1+£ {^). 



(129) 



For a fixed (3 < we have /3/y/rj —> — oo and we use the asymptotic expansion of D p (z) as 
z — > — oo which applies as z — > — oo (arg(z) = it) in (71), even if p is close to a positive integer. 
With p = 1 + e we have 



^i+ e (- 



2vr 



-e-2 



r(-i-e 

,1+s 



'-e z2 ' 2 [l + 0{z- 2 )] 



(130) 
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For e small we furthermore approximate T(— 1 — e) by e 1 + 0(1), which follows from the Laurent 
expansion of T(z) near z = —1. We thus obtain from (130) 

e^JW^) ~ ^= + ^v^Fee^ (131) 
and since term by term differentiation is permissible (see the discussion in Section 5) 



~ 1 + v^ee^^jp (132) 



Using (131) and (132) in (129) we see that e must be roughly of the order exp(— /3 2 /2r]); more 
precisely, 



-A 



E (ri)~--JL=e-*i[l-l3K{l3)] (133) 

and then —6 = r](l + e) leads to Proposition 4. 

We have thus shown that for f3 < there is a root of V = where 9 = 0{rf) with ~ —77. Now 
we show that there can be no roots in the range = o(rj), other than 9 = (which is an exact root 
for all 77 and /3). For j3 < 0, /3/y/rj — >■ —00 and we use (72) and set = 9*r] M . Thus 



D-,.,«- 1 (-|/3|/ v ^) ~ »p (^) rflf(^M=r) 

~M«-^f «,(£), < 134 > 

since 1/T(p) = p + 0(p 2 ) by Taylor series. Then using an analogous formula for D'_q, (■) and 
approximating D_q^m (— /3) ~ e~^ 2 ^ 4 and m{—@) ~ ^/Se - ' 32 / 4 we obtain 

V = e-PVe^-Wexp (^) v^F[l + o(l)]. (135) 

This asymptotic relation holds for 0^, fixed, 77 — > 0, and the o(l) error term holds uniformly on 
finite 9* intervals. But since the leading term is proportional to 0* we again conclude that 0* = 0, 
contradicting our assumption that a root has 9 = @(rj M ) for M > 1. If M = 1 the first asymptotic 
relation in (134) still holds and then we would find that V is proportional to 1/T(S) (with now 
9 = Sr]), which has roots at S = 0, —1, —2, .... The first root regains 9 = 0, the second is the 
spectral gap we analyzed above, and the roots 9 ~ —Nr] for N > 2 correspond to the higher 
eigenvalues. 

The expansion in (135) relies only on 9/rj being small. For any 9 = 0(77) we again obtain (135), 
with 9*rj M replaced by 9, and this excludes roots where 9/rj = o(l), except for 9 = 0. 



6.7 Proofs of Propositions 6 and 9 

We analyze V = for 77 — > and < (3 < /3*. We can no longer use (117) to approximate the 
parabolic cylinder function, as we will have 9 ~ — /3 2 /4. This corresponds to two saddle points in 
(114) (at i(/3± y/jPVM)) coalescing, see the discussion in Section 5 below (86). Now we must 
approximate D_q/ ^(fi I \Zv) an( i its derivative by Airy functions. We use the following proposition, 
which follows from [1, p. 689], and was discussed in (79)-(86). 
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Proposition 14. If A, B -too with A = -\B 2 + (\B) 2 / 3 5 and 5 = O(l), 

-B 2 /s/ 2 \A 



D_ A (B) = e 



B) 



4^ 1 \2' 



2\> 



^(1) ^(1) 



Ai(<y) + 



1 

*\2/3 

1 



(<5 2 Ai'(o~) - 2tfAi(<5)) + 0(5 



-4/3 > 



24/3^2/3 

T/ie error terms in (136) and (137) are uniform on finite 5 intervals. 
We let A = 6/n, B = /3/y/rj, 6 = x, and note that 



(136) 



(137) 



X 



rj An 22/3^1/3 



/3\2/3 



2/3, 



We rewrite V = as 



■iy-ei-fiD-g/niP/y/ri) = ^}D>_ eh (p/^j)D_ e {-f3) 



(138) 



(139) 



and recall that, by definition, /3* is the minimal root of D^ 2 0) = 0. For —9 = /3 2 /4 + 0(n 2 / 3 ) 
we use (136) and (137) in (139) and cancel some common factors to obtain 

■JW-0) + °(" 2/3 )l f Ai (x) + 2-" 3 /r 2/ y /3 (x 2 Ai(x) - 2Ai'( x )) + Of,, 2 ' 3 ) 



/3\V3 _ 



1/6 



5 /32/4 (-/3) + 0(r / 2 / 3 ) 



Ai'( x ) + 2- 4 / 3 /3" 2/3 " 1/3 (x 2 Ai / (x) - 2 x Ai( x )) + 0( V 2 / 3 ) 



0(r? 1/3 ). 



(140) 



The error terms are uniform on finite x intervals. The equation (140) applies both for < j3 < /3* 
and (3 ~ /3*, but its solution is different for these two cases. For < /3 < /3* the first factor in the 
left-hand side of (140) is 0(1), while it is o(l) if = & (or /3 « #,). 

First we consider < /3 < /3*. The right-hand side of (140) is O^ 1 / 3 ) so that x must be such 
that the left-hand side vanishes (to leading order in rj), which implies that Ai(x) = 0. Thus x must 
be close to a root of the Airy function, and the maximal root occurs at ao = —2.33810.... To 
obtain a more precise estimate we let x~ a o = r ? 1 ^ 3 Xi(/5, v) so that Ai(x) ~ f] l ^Xi{Pi r ?)Ai'(ao) as 
n — > 0. Then (140) becomes 



^ 2/4 (-/^ 1/3 [xiOM)Ai'(ao) - 2- 1 / 3 /3- 2 / 3 a Ai / (a ) + 0{r, 1 ' 3 ) 



/ 3 ^ 1/3 ?? l/3 



Ai / (a ) + O(n 1 / 3 )j Dp/^-P). 
Dividing (141) by n 1 ^ 3 and letting ij->0we conclude that Xi(fiiV) ~~ ^ Xi(P) as V — > 0, with 



-(«-^- 1/3 -(f)' /3 ^- 



(141) 



(142) 
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This leads to (18) and completes the proof of Proposition 6. 

We next consider (3 ~ /3*. Then Dpi 7^ and by Taylor series we have 

D' R2 , A (-P) 

' ( R , = ~m -&) + 0((J3 - &) 2 ), (143) 

where L is given by (25). Thus if we scale /3 — /3* = the left and right-hand sides of (140) 

are both 0{r] 1 ^) and we obtain the limiting equation 

-L-^-Ai( x ) = (/3,/2) 1 / 3 Ai'( x ), (144) 

which defines x implicitly in terms of W, and leads to Proposition 9. 

Finally, we show that when j3 > the equation V = can have no roots in the range 9 = o(l). 
We consider scales of the form 9 = 6*7] M with M > and exclude the possibility of roots that have 
= 0(rj M ). We write 



^D^-P)-^D^-P)^^ (145) 



and use (89) which for 9,rj = o(l) yields 



D _ e/ ^- ^ Wv ^ M (1 46) 
Also, using (67) gives 

£;D-o(-P) + ^D_ e (-(3) = 9D_ e -i{-P) ~ OD.^-fi). (147) 
dp 2 

Using (146) and (147) in (145) yields, to leading order in r], for 9 = 9*r] M , 

V = e^ M D_ e ^ (^) [/T + D_i(/3)] [1 + o(l)]. (148) 

For M > 1 and r] — > we can approximate D_q^m-i{^=) ~ exp[— (3 2 /(4i])], for M = 1 we have 
D_Q tT] M-i{-^) ~ exp[— (3 2 / (4r])](^/rj / (3) e * , while for < M < 1 we must approximate the parabolic 
cylinder function using the result in (117), which applies for large index and large argument. But in 
all cases the approximation leads to V in (148) being proportional to 9* with a positive multiplier. 
Thus we again conclude that 9* = 0, contradicting the existence of the root(s) where 9 = Q(r] M ) 
for any M > 0. Note that unlike /3 < 0, the scale 9 = 0(ry) does not lead to roots, but only a 
change in the expansion of D_Q itr} M-i{-^=). 

In obtaining (148) we used only the fact that 9 = o(l). For any 9 = o(l), (148) holds, with 
9*r] M replaced by 9. Again, the expansion of D_qj ti {P/ y/rj) will be different as 9/rj — > 0, 9/rj — > 00, 
or 9/r] = 0(1), but the first multiplicative factor in (148) (= 9*rj M = 9) shows that there can be 
no roots in any range where 9 = o(l), except for the root at 9 = 0. 
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7 Monotonicity of the spectral gap 



The surface sketched in Figure 1 suggested certain monotonicity properties of r(/3, 77), and these 
were partially confirmed by the various asymptotic formulas in Section 2. We now establish these 
analytically, for all values of (/3,rj). We shall obtain: 

Proposition 15. Let sgn(z) = +1 if z > 0, sgn(z) = —1 if z < and sgn(0) = 0. Then 



( Tjj^j = -sgn(?7 



1) 



(149) 



Hence, for a fixed rj < 1 the spectral gap r is an increasing function of (3, while it decreases with f3 
for fixed rj > 1. If rj = 1, r(/3, 1) = 1 is constant. 

To establish this result it is useful to set 



V(P;r ] ,f3) = V(9;r ] ,f3), 9 = -P 



and then in view of Theorem 1 

V(r(P,7i);Ti,P) 
By implicit differentiation of (151) we have 



0. 



dV 



dr 
dp 



dV 
~dp 



0. 



(150) 



(151) 



(152) 



By definition, r is the minimal positive solution of V(P; rj, (3) =0 and we also note that V(0; rj, (3) = 
0, as 9 = is a simple pole of p(x; 6) in Theorem 10, which corresponds to the steady state limit 
in (6). Thus P = and P = r are consecutive zeros of V(P;rj,f3) = 0. To determine the sign of 
dr/d/3 in (152) requires that we know the signs of dV/dP and dV/d(3 at P = r. For the former we 
can compute dV/dP using the expression in Theorem 1 and the integral respresentations in (61), 
(64) and (65). However, an indirect argument leads immediately to the value of sgm{dV/dP\ P=r ). 

The solutions of V = for P > correspond to poles of the Laplace transform p{x\ 9) in 
Theorem 10 and these are the eigenvalues Aw for N > 0, with Ao = and Ai = r. From the 
general theory of the one-dimensional Schrodinger equation, the eigenvalues are all simple (see the 
discussion below (52)) and the equation V = has simple zeros. Hence, dV/dP\ p=r ^ 0. We also 
note that if V had, say, a double zero at P = r, then (28) would imply that the spectral expansion 
of p(x,t) would involve the terms e~ Xlt = e~ rt and also te~ rt , and this would contradict the self- 
adjointness of the Schrodinger equation in (49). Now, since P = and P = r are consecutive 
simple zeros on the real axis of the entire function V (as a function of P) we must have 

dV \ fdV 
= — sgn I 



sgn 



dP 



V dP 



p=o 



(153) 



Computing the right-hand side of (153) is much easier than computing the left-hand side, as we 
show below. 



We define the functions I(z) and J(z) by 



P=0 



Br 



{\ogu)e- zu e u2/2 du, 



(154) 



J(z) 



, u2 / 2 du. 



Br 
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Then we expand V in (150) in Taylor series about P = to obtain 
V=-y/fj [D (-P) + PI{-0) + 0(P 2 )] 



P + 0(P 2 ) 



Dr 



+ i 



+ 0{P 2 



Using Dq{z) = e z2 / A : D' (z) = — 7j e we obtain from (155) 



dV 



cxp 



P=0 



-/3 2 /4 



j 



But an integration by parts shows that 



zl{z) - J(z) 



i\^2K J Br 

e 



(z - u) log(u)e- zlt e u2/2 
y log(u) d (- e - 2U e u2/2 ) 



i\/2TT J Br u 



°?l*du 



-e/2 



(155) 



(156) 



(157) 



Here we also used the fact that a parabolic cylinder function of order P = —1 can be expressed in 
terms of the standard error function, or probability integral. Using (157) in (156) we have 



dV 



p=o 



exp 



4 47/ 



-e/2 



d£ + exp 



/3 2 /3 2 
+ 



4 4r ?/ hi^\ 



< 2 ' 2 di 



(158) 



exp 



4t/ 4 



d£ + 



K e -ei2 di 



e 



so that dV/dP\ p=0 < and hence, in view of (152) and (153), 



sgn 



dr 

df3 



sgn 



dV 
~d/3 



P=r 



(159) 
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Now, 



dV{P\y,P) 

d/3 



(160) 



D P {-P)D'^ h 



-f Ml 



-P)D P/V 



Using the parabolic cylinder equation D p (z) 



P — t;)Dp(z) we can simplify (160) to 



dV 



K] — 1 
V 



p/v 



JL 



+ 



P 2 



P D P {-P)Dp h 



p_ 



• (161) 



When P = r we can further use the fact that V(r; rj, /3) =0 to simplify the right side of (161). We 
will need to separately consider the two cases D r {—(3) = (a degenerate case that occurs rarely) 
and D r (—f3) ^ (which is typical). 
In the degenerate case we have 



sgn 



dV 

d/3 



P=r 



sgn(r / - 1) • sgn {D' r (-/3)) • sgn \ D' r / n (q=)) > ( 162 ) 



and thus 



sgn ( t^J = - sgn(?7 - 1) sgn (D' r (-(3)) sgn ( 



(163) 



We note that if D r (—f3) = then certainly D' r {—j3) / 0, as discussed in Section 5 below (69). But 
if both V = and D r (—f3) = then certainly D r / v (f3/y/rj) = 0. In the non-degenerate case we can 
rewrite V = as 



r/r/ 



D r / V I — — 



P_ 



(164) 



Using (164) in (161) and (159) we conclude that 



sgn 



/ dr 



sgn(ry - 1) • sgn(£> r (-/3)) 
•sgn (r> rh ^g 11 



(165) 



\ Dr (-P)f + (D' r (-P)) 2 



We proceed to determine the signs of the various terms in (163) and (165). 

It proves useful to understand the behaviors of r(/3, rj) as j3 — > =too for a fixed rj. By a calculation 
completely analogous to that used to establish Proposition 4, we find that 



to \ (1-V)y/V 
r{p, rj) — rj ~ - — / 1_ exp 



-oo. 



(166) 
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Whereas Proposition 4 applies for rj — > with fixed P < 0, (166) applies for fixed 77 as P — > — oo. 
By expanding (14) for /3 — > — oo and (166) as 77 — > we see that the two agree in this intermediate 
limit. Thus for (3 large and negative, r is exponentially close to 77, as could be expected since then 
almost all of the probability mass in the model migrates to the range x > 0. A similar analysis as 
P — > +00 shows that 

77 — 1 



r(P,v) ~ 1 



-/? 2 /2. 



/3V27T 



7/^1, /3^+oo, 



(167) 



which can also be obtained simply by using (166) and the symmetry relation in (13). For /3 large 
and positive the probability mass migrates to the region x < 0. Note that (166) and (167) suggest 
that dr/d/3 has the oppositive sign as 77 — 1, at least for \(3\ sufficiently large, and this we proceed 
to establish for any f3. 

Returning to (165) we proceed to show that D r (—[3) and D r / 7] (/3/^) always have opposite 
signs. In general, suppose that we have two real analytic functions F{x) and G{x). The ratio 
F(x)/G{x) can only change sign at a point x* where F(x*) = with G(x*) 7^ 0, or G(x*) = 
with F(x*) 7^ 0, or, possibly, where F and G both vanish but have zeros of different orders. Thus 
if F(x) and G(x) are non-zero, or if their only zero(s) coincide and they are of the same order, 
then F(x)/G(x) cannot change sign. Then determining sgn(F(x)/G(x)) requires only that we 
evaluate the ratio at a particular point, which could be x = ±00. But, we showed below (163) 
that D r (—f3) and D T i rj {f3/ \/rj), as functions of (3 for a fixed 77 > 0, can only vanish simultaneously. 
Furthermore, if D r {—(3) vanishes at some /3 C , then for f3 near j3 c but (3 7^ (3 C we can rewrite the 
equation V(r; 77, (3) = as 

D r (-P) D' r (-f3) 



Then letting /3 — > {3 C leads to 



lim 



D r (-P) 



D rh {-(3/^j) 



(168) 



(169) 



where r c = r(/3 c , 77). But (169) shows that D r (—/3) and D r / TI (P/^/rj) must vanish to the same order 
at /3 = f3 c (in fact they must have simple zeros there). Thus we conclude that D r (—(3)/D r / v ((3/^/rj) 
cannot change sign. To determine this constant sign we can let either (3 —> +00 or /3 — > —00 as 
then we have asymptotic formulas for r. Using (167), (72) and the fact that Di(—/3) = —f3e~@ 7/4 
we find that 

Dr/nW/Vv) 1-1 1 

—p 7 ! rj 2r > exp 



D r {-(3) 



-)/3 2 
77 



+00 



(170) 



and thus D r / v (/3/y/rj) and D r (—f3) have opposite signs as /3 — > +00, and thus this is true for all f3. 
Note also that for (3 —> —00, (166) and (70) lead to 



Dr/nW/y/r}) 
D r (-P) 



■(-/?) 



1-7?. 



1 



exp 



1 



1 



1 



P 7 



P 



-00 



(171) 



and this verifies the conclusion about opposite signs. We have thus simplified (165) to 



Qt \ 

sgn ( tt^ J = - sgn(?7 - 1) sgn 



^)W-/3)) 2 + (^(-/3)) 2 } (172) 
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in the non-degenerate case. In the degenerate case we conclude from (169) that D' rc {—j3 c ) and 
D' r ^(Pc/y/v) have the same sign, and then (163) shows that sgn((9r / 'd/3) = — sgn(r/ — 1), which 
establishes Proposition 15. 

It remains to show that the last factor in (172) has always positive sign. Let us define 

U(P,z) ={P-^) Dl(z) + [D' P (z)]\ (173) 

and we consider H as a function of both P and z. We clearly have H(P, 0) = PDp(0) + [D' p (0)] 2 > 
for P > 0, with H(0,0) = 0. Also, H(Q,z) = for all z, in view of (63). We consider P > and 
z > 0. 

We shall show that T~L(P,z) > for all z > when P > 0. For z — > oo the estimate in (89) 
leads to 

(D'p(z)) 2 z 2 - P 

and then from (173) and (70) 



H(P, z) ~ ~ Pz 2p - 2 e" z2 / 2 , z -> +oo (175) 



z 2 



so that "H is positive for z sufficiently large. By differentiating (173) with respect to z we obtain 

^ = 2D P {z)D' P (z) (j° " - 1 + 2D' P (z)D' P (z) (176) 

= - |iP 2 (z) - D P {z)D' P (z) 
= - PDp{z)D P ^{z). 

Here we also used (60) and the recurrence (67). From (176) we conclude that % has maximum 
or minimum values at roots of Dp{z) and Dp-i(z), as functions of z. As discussed in Section 5, 
there are at most finitely many of these. But if Dp(z*) = for some z* then D' p {z*) ^ and 
n{P,z*) = {D' p (z*)} 2 > 0. If D P _i{7) = for some z then (173) and (67) show that 

H{P,z) = (P~ Z -P) D 2 P & + [D' P &] 2 (177) 
= PD P (z), 

which is again positive for P > 0. Note that we cannot have simultaneously Dp_i(z) = = Dp(z), 
for then (67) would imply that D' p {z) = also. We have thus shown that ~H(P,z) is (for P > 0) 
positive at z = and as z — > +oo, and also % > at any maximum/minimum value of T~L. We 
then conclude that 

n(P, z)>0 for P > and z > 0. (178) 

Note that if % becomes negative at some z = z 1 then % would need to reach a minimum value at 
a point z" where T~L < 0, since for sufficiently large z we again have % > 0. 

Now we let P = r(/3, rj) > and z = —(3 and use (173) and (178) in (172) to conclude that 

sgn f |0 = - sgn(? ? - 1), < 
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and we have thus established Proposition 15 for /3 < and all i] > 0. To show the result holds also 
for (3 > 0, we need only use the symmetry relation in (13), which shows that if r increases with 
/3 for f3 < and < rj < 1, (resp. r/ > 1) then r will decrease with (3 for /3 > and r/ > 1 (resp. 
< r? < 1). 

Alternately, we can use the relation (168) (in the non-degenerate case) and (172) to conclude 
that 



sgn 



dr \ 
dp) 



sgn(rj — 1) sgn 




D r/V 



+ r] 



(179) 



and apply (178) with P = r/r] > and z = (3/ ' ^ffj > 0. This concludes the proof of Proposition 15, 
which was suggested by our numerical and asymptotic results. 
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A Appendix A 

We now prove Proposition 8. First consider 7 — > —00. Since (20) is equivalent to -D#_i(7) = 
for R 7^ 0, the asymptotic formula in (72) shows that Dr_\{^) is positive for 7 — > —00, unless 
r(l — R) is singular, and the minimal singularity occurs at R = 1. For R close to 1 we must use 
the asymptotic formula in (71), with arg(z) = arg(7) = tt. Then Dr_i(^) = implies that 

e -7 2 /4 7 fl-i „ ^(_ 7) -V 2 /4 7 -> -00. (180) 

1 (1 - K) 

Then using T(l - R) = (1 - R)- 1 + 0(1) as R -> 1, (180) leads to 

_ 7e -7 2 /2 

R - 1 ' , 7 -> -00 (181) 

V27T 



which is the result in (21). 

Now consider the limit 7 — > +00. For a fixed -R the asymptotic formula in (70) shows that 
Dr-i^) ~ 7 iJ ~ 1 e -7 / 4 , 7 — > +00, which is strictly positive. Thus to capture the zeros of Dr_i(7) 
in this limit R must be itself large, so that we enter the oscillatory range of the special function. 
As discussed in Section 5 the transition to oscillatory behavior occurs when R ~ 7 2 /4 and then 
we can approximate Dr-i^) by Airy functions, with the leading term given in (86). Thus with 
R = 7 2 /4 — (7/2) 2 / 3 <5 and 7 — > 00 the minimal root corresponds to the maximal root of Ai(J) = 0, 
which occurs at 5 = clq, leading to (22). 
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B Appendix B 



We discuss the singularities of (28)-(30) in the complex #-plane and thus establish Proposition 11. 
As discussed in Section 5, D-g(f3) is an entire function of 9, so that the only singularities of (28) 
are the zeros of V{9; r], 0). The existence of an infinite sequence of zeros and the fact that they lie 
on the real axis (Im(#) = 0) follows from standard ODE theory, which was discussed in Section 4. 

Now consider (29) and (30). The factor V{9) has simple poles at 9 = 0,-1,-2,-3,.... If 
9 = —M, M > we can simplify V by using (63), so that 



V(-M; r,,())=- ^Re M (-^ 2 / 4 D' M/rj (JL^J (182) 



+ dp L 



Dm/t) 



Similary, (27) and (63) lead to 

M(—M; rj, 0) = ^Re M ((3)e-^^D' M/v (JL^J (183) 

He M (-/3)e-^/4l 



_d_ 

dp 



Then we use HeM( — 0) = (— l) M HeA/(/3), as the Hermite polynomials are odd/even functions 
according as M is odd/even, and comparing (182) to (183) we find that 

V(-M; v , (3) = {-l) M+l M(-M; V , (3); M = 0,1,2, ... . (184) 

From (184) we conclude that either V and M are both zero at 6 = —M, or neither is zero. If 
9 = -M and V(-M;r),0) ± 0, then M/V = (-1) M+1 . But then 

Dm(xo+ n + Dm{ . X0 

= e-<' J+ «> ! / 4 [He M (i„ +P) + (-l)" +1 He„(-i -/?)]= 

so that the last factor in (29) vanishes, and thus 9 = —M is not a pole of (29). Similarly, 9 = —M 
will not be a pole of (30). This shows that 9 = —M can only be a pole of (29) and (30) if 9 = —M 
and V and M simultaneously vanish. 

Conversely, suppose that V and Ai both vanish, say at some 9 = 9*. First we assume that none 
of the three D-$ t (— 0), D-q 9 {0) and D_g t / V (f3/ ^frf) are zero. Then the equations M. = V = may 
be rearranged to give 

#_,.(-/?) _ 2? fl . Q8) _ 

D_ e ,{-(3) D_eM VV D^ /r) ((3/ ^j)' 1 ) 

But the first equality in (185), along with the Wronkskian identity in (68), leads to 



= -D'_ e S-{i)D-oM ~ D^{-P)D'_ e ^) = ^r. (186) 
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But then r(#*) must be infinite, which leads us back to the case —9 = M = 0, 1,2. . . which we 
already discussed. 

Finally suppose that D^g t (—/3) = 0. Then certainly D'_ e (—(3) ^ 0, and V = implies that 
D-9„/n(l3/y/v) = 0) an d then Ai = leads to D_Q t (f3) = 0. But then (186) again leads to the 
conclusion that 9* = —M = 0,-1,-2,.... Starting with the assumption that D_e t (/3) = or 
D-o^/^iP/y/rj) = leads ultimately to the conclusion that all three denominators in (185) must 
vanish, and then again (186) leads to 9* = —M. 

We have thus shown that simultaneous zeros of V and Ai can occur only if 9 = 0, —1, —2, 
Indeed this does occur precisely when 9 = 0, which corresponds to the steady state limit p(x, oo). 
In all cases we showed that a singularity of (29) or (30) necessarily has V = 0. Thus the equation 
in Theorem 1 captures all of the singularities of p(x; 9). 



C Appendix C 



We consider the expressions in (28)-(30), in the limit ?] — > + , and thus establish Proposition 
13. In this limit we can simplify D_ e / ri ((r]x + /3)/y / ?7) using the asymptotic formula in (77). We 
can also use (64) to obtain an analogous formula for D'_q,((3/ y/rj). Let us assume first that 9 
is positive and real. The expansion (77) follows from a saddle point approximation to (61), as 
discussed in (77) and below. The expansion of D' p (z) for z,p — > oo with z 2 /p > 4 is the same as 
that of D p {z), except that the factor z/2 — u in the integrand in (64) becomes frozen at the saddle 

4p] . It follows that 



u 



zv + 



D p (z) 



4p; z,p — >■ oo, — > 4 
p 



(187) 



and hence (setting z = (3/ 'y/rj and p = —O/rj with 77 — >• 

(^) 



D_ 



A similar argument shows that 



eh 



(3 + rjx 



-V^+40 as T] -+ H 



(188) 



as 7] 



0+. 



(189) 



Then writing V as 



V 



-D. 



■oh 



JL 



D'_ e {-P) + ^ 



D_ 



-e/v 
eh 



(190) 



D_ 



eh 



-D'_ e (-P) + -v / /? 2 + 4^_ 9 (-/3) 
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we see that as r/ — > the expression in (28) becomes that in (39). We also have 



M D'_ e {P) + ^6 + ^/±D- e {f]) 

— > : , 77 — > (191) 

V -^(-/9) + ^ + /3V4^(-/3) 

which can be used to obtain the limits of (29) and (30), and this agrees with the results we 
obtained in [23]. Throughout this calculation we divided several times by D_g/ 71 (/3/^/rj), which 
was permissible since, for 9 > 0, we are outside of the oscillatory range of the special function, as 
we discussed in Section 5. 



D Appendix D 

Here we discuss the discrete M/M/m + M model. We shall obtain an explicit, albeit complicated, 
expression for the Laplace transform (over time) of p n (t) = Prob[iV(t) = n \ N(0) = no], where N(t) 
is the number of customers in the system. Then we will give an alternate derivation of Theorem 1, 
by evaluating the discrete model in the limit m — > oo with p = A/ u = m + 0{^/m). The analysis 
here closely parallels the proof of Theorem 10, so we just give the main points. 

We solve the following infinite system of ODEs (we assume time is scaled to make the service 
rate p, = 1, so that p = A): 

p'S) = Pi(t) - PPo(t) (192) 

Pn(t) = P [Pn-l(t) - Pn{t)\ + (n + lK+l(t) ~ np n (t), 1 < 77 < 777 - 1, (193) 
p' m (t) = p\p m -l(t) -Pm(t)\ + (777 + T])p m+1 (t) - mp m (t), (194) 



Pn(t) = P \Pn-l(t) ~ Pn(t)} + [m + (n - 777 + 1)7/] Pn+l(t) (195) 
— [777 + (77 — 777)77] p n (t), 77 > 777 + 1 

with the initial condition p n (0) = 5(n,no). We need to consider the cases uq < m and uq > m 
separately, as the discrete model has no analog of the symmetry relation in Proposition 3. 

Introducing the Laplace transform p n {9) = / °° e~ et p n {t) dt, we first consider the case < no < 
777 and then (192)-(195) become 

= Pi-(p + 0)p (196) 

- S(n,n ) = ppn-i - {p + 6 + n)p n + (n + l)p n +i, 1 < n < m - 1 (197) 

= pp n -i - [p + 9 + m+ (n - m)rj\ p n + [?n + (n - in + l)n] p n +i, n > m. (198) 

We solve (196)-(198) using a discrete Green's function approach. We begin by introducing the 
functions F n , G n , H n , I n ; these are defined by the contour integrals 

p n ~ l (9 + £-l)(9 + £-2)...(9 + l)8 
~ ^ (n-e)\ I\ ' 
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Figure 5: A sketch of the branch cuts and the contours C\ and 



H n {6) 



1 

2vr7 
1 



C 



z~ m/r > dz 
Z - mh dz. 



(201) 
(202) 



m / r? exp [ — imrj 1 arg z\ 
arg(z — 1) — 7r]} (so that 



Z^O/v ^n+l—m 

Here Cq is a small loop about z = 0, C\ goes from —00 — is to —00 + ie, encircling z = 1, while C2 
goes from — 00 — to — 00 + ie encircling z = 0. The contours Ci and C2 are sketched in Figure 5, 
and in (200) we use the branch (z — l)~ e = \z — l\~ e exp[— 16 arg(z — 1)] with — ir < arg(z — 1) < vr; 
in (201) (z-iy 6 ^ = \z-l\- e ^exp [- iOrj" 1 arg(z - 1)} and z~ m ^ = \z\ 
with -vr < argz < vr; and in (202) (1 - z)~ e ^ = |1 - zj-^exp { - iQif 1 

l-e/n = 1 for > 0). Note that if 77 = 1 then G n = H n and F n = I n (then the contour C 2 in (202) 
may be deformed to Co). 

We can easily verify that F n and G n satisfy the homogeneous form of (197) (with <5(n, no) 
replaced by zero) and thus give two linearly independent solutions of this difference equation. 
Similarly, H n and J n give two solutions to (198). We use the functions in (199)-(202) to construct 
p n (6), making use of (194) (or (198) with n = m) to relate the ranges n < m and n > m. The 
analysis is completely analogous to the proof of Theorem 10 in Section 6, so we give below only 
the final result 
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Here we suppressed the dependence of F n , G n , H n on 6. 
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The above holds for all < no < m, and if uq = m (starting with all servers occupied but an 
empty queue), (204) is not needed, and then p n (6) somewhat simplifies, to 



Pn{0) 



F m H n , n>m 



F m H m -\ — H m F m _i I H m F n , < n < 



m. 



(206) 



Note that (203), (204) and (205) are similar in form to (28), (30) and (29), respectively, with 
H m G m -i — G m H m _i and F m H m -i — H m F m -i playing the roles of M. and V. The functions 
in (199)-(202) are all entire functions of 0, and the singularities of p n (Q) are determined by the 
equation 

F m (e)H m ^(e) - H m (e)F m ^(e) = o. (207) 

For initial conditions no > m we need to solve (196)-(198), but now with the left-hand side of 
(197) replaced by zero and the left-hand side of (198) replaced by — <5(n, no). For uq > m the final 
result is 



Pn(0) 
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(208) 
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(209) 



Ue) = -( p - 
p \v 



m—riQ 



r(no + 1 — m + m/rj) 



H F 

11 no ± n 



T(l + m/rj) F m H m _i - H m F m _i ' 
< n < m. 



(210) 



When no = m, (209) is not needed and then (208) and (210) reduce to (206). Again the singularities 
are determined by (207). 

Now we evaluate these results in the limit of m — > oo with the scaling 



n = m + \fmx, no = m + \/mxo, p = m 



m/3, 



(211) 



where x, xq, f3 are 0(1). We shall thus give an alternate derivation of Theorems 1 and 8 of Section 2. 
Let us also scale z = 1 — £,/ y/rn in the integrands in (199)-(201). Noting that 

pz — n log z = (m — y/mf3) ( 1 ^= J — (m + \frnx) log (1 ^= ] 

V V m J V v m J 
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we obtain a limiting form of (199): 



F n {9) ~ ^=^=. [ rV*+^ 2 / 2 d£ (212) 



m e»/2 e p 



\/2irm 



e-^ + Wl*D_ e (-x-P). 



Here we used the integral representation in (61) for the parabolic cylinder function D. In (212) 
Br is a vertical contour with Re(£) > 0, which can be used to approximate C$ in (199) with this 
scaling of z. A completely analogous expansion of (200) leads to 



m e/2 eP 



G n (9) ~ '-^—Le-^r/* D _ e{x + j8)j (213) 
V27rm 



and, after some calculation, we obtain from (201) 



Using (201) and Stirling's formula we also have 



(??x + f3) 2 
4r] 



D-e h ( ^ ) ■ (214) 



V^. p m- no-1 „ l e ^o e ^/2_ (215) 

m! m 

Next we consider the limiting form of (207). Noting that F m — F m _i can be computed from (199) 
by multiplying the integrand by 1 — z and setting n = m, we have 

F m H m -i — F m -\H m (216) 

— H m (Frn Fm—l) F m (H m H m —\) 

V zirm 

e 

„ J_ m */2 ( ™\ 2 " e P(l+l/.) e -/3 2 /4 e -/3V(4r ? )^ V(0;r?)/3) 

where V is as in Theorem 1. Here we also used (212), (214) and (68). 

An analogous approximation to H m G m -i — G m H rn _i can be used to show that 

H m G m -i - G m H m -\ ^ M(9; rj, P) 

F m H m -l - H m F m -l V(9] 7?, 0) ' 



(217) 



Then using (212)— (215) and (217), we see that the expressions in (203)-(205) reduce to those in 
(28)-(30), up to a factor of l/^/m in the former, which arises due to the fact that p n (t) is normalized 
by a sum over n while p(x, t) is normalized by an integral over x. We have thus given an alternate 
derivation of Theorem 10. Note that all of the asymptotic calculations do not involve scaling time t 
or the transform variable 9. 
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Finally, we discuss the uniformity of the approximation in (216) for small values of 9. This will 
show that A(0) = A(6;m, f3, 77) = F m H m -\ — F m -\H m can have no roots for m — > 00 in ranges 
where 9 = o(l). We set p = m — fi^/m and consider finite intervals of /3 and 77, with 77 > 0. First 
we note that, using (199) and (201) (with n = m and then n = m — 1) 



F m (0) = — v H m (Q) 
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T(rn/rj + 1) 
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(218) 



It follows that F m (0)i? m _i(0) = F m _i(0)iT m (0) and thus A(0;m,/3,r/) = for all values of m,(3,rj. 
The pole at 9 = corresponds to the steady state limit, which exists for all m and f3, for 77 > 0. We 
expand A for fixed finite j3 and fixed 77 > 0, as m — > 00, which will refine the leading order result 
in (216) and show that the higher order terms remain smaller than the leading term, for 9 = o(l) 
as m —> 00. Setting z = 1 — £,/y/m in the integral in (199) leads to 



m 



~ e/2 F m (9) 
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(219) 



where C is the image of the contour Co and 
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(220) 



where fj(^) is a polynomial in f of degree 3j. We have Pi(£) = f + ±£ 3 , P 2 (0 = £ 2 + n£ 4 + t^ 6 , 
etc. Using (220) and some contour deformation (as \z\ < 1 in (199) implies asymptotically that 
Re(£) > 0), we obtain the asymptotic series 
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- e/2 F m (9) 



e-^D-ei-fl + J^m-WfMf}) 



(221) 



where 



fj(0,l3) 
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2ni 



Br 4 



my 



(222) 



and Re(£) > on the vertical contour Br + . Since the Pj are polynomials, the integral in (222) is a 
finite sum of parabolic cylinder functions of different orders (for example f\ involves D\_q{-) and 
Z?3_0 (- ) ) , and thus each fj is an entire function of 9. A completely analogous calculation shows 
that 
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Next consider H m (9) (setting n = m in (201)) in the same asymptotic limit. Now we scale 
obta 

V T] / 



z = 1 + £ yfr]/m and obtain 

> m \-0/(2r)) 
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e e/^ r o n{c . m)dt 



(224) 



with 



exp 
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mlo£ 
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where Pj are again polynomials in £. Then we obtain the asymptotic expansion of (224) as 
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(226) 



(227) 



so again hj is a finite sum of parabolic cylinder functions of different indices, all with argument 
—(3/^/rj. A completely analogous calculation shows that 
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(228) 



Combining (221)-(223) with (226)-(228) leads to a refinement of (216) into a full asymptotic ex- 
pansion of A(9;m, f3,rj) in powers of m -1 / 2 . Then we can replace V in (216) by 



where 



V(0;77,/3)+I>- J ' /2 V#;7/,/3) 

3=1 



V,(0; rj, fl) = Y^[h k (9, 0, v)fj-k(0 + h k (9 - r), /3, n)fj- k (6, /?)]. 

fc=0 



(229) 



(230) 



Here Vo = V and /o and ho can be identified from (221) and (226). For (230) we see that each 
Vj is a finite sum of products of two parabolic cylinder functions, possibly of different indices and 
arguments. Thus each Vj is an entire function of 9, and the expansion in (229) is uniform in finite 
fl,9,r] intervals, with 7] > 0. Setting 9 = and using the fact that A(0), which is itself and entire 
function of 9, vanishes for all values of m,(3,r], we conclude that Vj(0;rj,j3) = for each j. From 
(158), V has a simple zero at 9 = and then Vj for j > 1 have zeros of orders > 1. Thus for any 
9 = o(l) as m — > oo the correction terms in (229) remain smaller than the leading term, and thus 
A(#) cannot have a zero for 9 = o(l), except the one at 9 = 0. We have thus shown that for fixed 
j3, fixed r\ > 0, and to — > oo the zeros of A in (216) in the range 9 = 0(1) must approach the zeros 
of V, and that A has no zeros where 9 = o(l), except for 9 = 0. 
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E Appendix E 



Here we establish the positivity of the correction term in (15). We thus take ft > ft* and show that 
A(ft) > 0. 

We begin by noting that V(0, ft) = and V(r (ft), ft) = 0, and since r (ft) is the minimal positive 
solution of V(p, ft) = it follows that and ro(ft) are consecutive zeros of V(p,ft). These zeros are 
necessarily simple, due to the general results for the one-dimensional Schrbdinger equations that 
we discussed in Section 4 (see the discussion surrounding (53)). It follows that 



fdV 
Sgn W 



p=0 



fdV 
3gn W 



p=r (/3) 



(231) 



The right-hand side of (231) appears in A(ft) in Proposition 5, but the left-hand side is much easier 
to determine. Using (124) and the definition of V in Proposition 5, we have 



dV 







dp p=o dp 



^[D'^-ft) - D p (-/3)y/p/I=pl 
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(232) 
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A)(-^hvW ~P + D' p (-ft)/D p (-ft)} 
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+ e 
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e~ u2 / 2 du e"^ 4 > 0. 



Here we also used the facts that f3 > and V(0, j3) = 0. From (231) and (232) we conclude that 



fdV 
sgn — 
V op 



p=r o (0) 



(233) 



If we can show that D ro (-(3) < for all (3 > (3* then (17) shows that A(J3) > 0. 
Consider = D 



'r (S){~P) as a function of (3. As discussed in Section 6.3, T((3) is an 
infinitely smooth function of f3. This function cannot change sign, for if J-((3 C ) = for some f3 c 
then D p (-/3) = for p = r (/3 c ) and also, since V(r (P),P) = 0, D' p (-(3) = 0. But D p and D' p 
cannot simultaneously vanish, as discussed in Section 5 below equation (69). To determine the 
constant sign of J- (ft) we need only evaluate this at one particular point. For example, when 
ft = 2 then r$ ~ .9323 and D ro (p\(—ft) ~ —.8275. We can also let ft — > oo and use the facts that 
r (ft) -> 1 and D ro{(3) (-ft) ~ -fte~^^ < 0. Thus sgn(D ro (ft)) = -1, and (17) and (233) show 
that sgn(A(ft)) = +1. 
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